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The  CO  fundamental  vibration-rotation  spectra  resulting  from  the  interaction  of  discharged 
nitrogen  with  carbon  monoxide  at  low  pressure  ( ~3  mTorr)  exhibit  bimodal  rotational 
distributions.  We  have  identified  14  vibrational  levels  of  a  rotationally  relaxed  (80  K) 
component  and  eight  vibrational  levels  from  a  rotationally  excited  component.  The  eight 
rotationally  excited  bands  are  best  reproduced  by  a  statistical  distribution 
£  =  £*  +  Ev  *>  3.7  eV,  which  provides  sufficient  population  in  the  region  of  the  Fortrat 
reversal  (7~90)  to  account  for  the  observed  £ -branch  bandbead  formation.  The  rotationally 
relaxed  vibrational  levels  are  populated  by  single-  and  two-quantum  transfer  from  N2(u), 
N2(t>)  +  CO— N2(t>  —  1,2)  +  00(i>  =s  1,2),  and  radiative  cascade  from  CO(.4)  produced  by 
quenching  of  AT2(a')»  N2(n'  *2^ )  +  CO— N2(£,0)  +  CO  (A  *11)— CO(0<9)  +  hv,  and 
relaxation  of  the  rotationally  excited  component.  Kinetic  and  energetic  arguments  indicate 
that  a  branch  of  N2 (o')  quenching  N2(a'  '2  ~  )  +  CO— fi2(X,v)  +  CO(0<14,7)  is  responsible 
for  the  rotationally  excited  component.  Surprise!  analysis  indicates  two  dynamic  mechanisms 
are  responsible  for  the  rotationally  excited  component.  We  have  modeled  the  vibrational 
distribution  of  the  rotationally  excited  component  with  equal  contributions  from  a  statistical 
(all  v)  process  and  a  process  favoring  excitation  of  low  vibrational  levels  (u<4). 


INTRODUCTION 

Modest  energy  depositions  (typically  <  1  eV)  in  CO 
interna]  states  have  been  observed  in  photochemical  and  ab¬ 
straction  reactions.  In  general,  large  fractions  of  the  energy 
above  threshold  are  manifested  in  product  internal  states. 
For  example,  photolysis  of  OCS  at  157  nm,1  acetone2  at  193 
nm,  and  H^CO3,4  produces  vibrationlly  excited  CO  with  ro¬ 
tational  excitations  up  to  0.9  eV.  Studies  of  energy  partition¬ 
ing  in  CO  from  the  reaction  of  hot  H  atoms  with  CX>2  indi¬ 
cates  nearly  1.0  eV  in  CO  internal  states’  with  the  rotational 
distribution  in  the  v  =  0  level  following  a  statistical  model. 

Greater  energy  depositions  are  observed  for  energy 
transfer  reactions  since  these  interactions  have  little  or  no 
threshold  energy.  Quenching  by  0(2  '/)),*  Na(3  2P),7* 
1(5  ^./s),’  and  Br(42f>I/2),v  and  Hg(63/»,  and  63P0)'° 
have  been  observed  to  produce  vibrationally  excited  CO.  Ro¬ 
tational  excitations  up  to  0.8  eV  have  been  reported  for  the 
Na(3  2P)  +  CO  quenching  reaction,8  which  accounts  for  a 
large  fraction  of  the  total  2.1  eV  exoergicity. 

The  excitation  process  responsible  for  the  CO  excitation 
reported  in  this  paper  is  energy  transfer  from  a  metastable 
state  of  nitrogen.  Energy  transfer  to  CO  from  discharged 
nitrogen  has  been  examined  closely,  particularly  with  re¬ 
spect  to  the  N2-CO  laser."12  High  CO  vibrational  excita¬ 
tion  is  produced  from  near-resonant  N2(u)  energy  transfer 
and  subsequent  CO(u)  collisional  up  pumping13,14  (reac¬ 
tions  (la)  and  (lb)] 

N2(u)  -f  CO(u')— N2(0—  1)  +CO(0'  +  1),  (la) 

CO(0)  +  CO(0')-CO(0-1)+CO(0'  +  1).  (lb) 
Although  high  CO  vibrational  levels  are  produced  in  this 
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manner,  no  rotational  excitation  is  observed,  which  is  con¬ 
sistent  with  the  small  energy  defect  for  the  near-resonant 
process. 

The  energy  defect  from  N2(n  'll )  quenching13 

n2(o 'm  +  co-n2(a»  +  cou'n,i/,;')  (2) 

has  been  determined  to  be  largely  manifested  in  rotation,  but 
the  total  energy  defect  is  small  (  <  1000  cm1),  much  less 
than  the  rotational  excitations  observed  here.  The  quench¬ 
ing  of  N2(4  32„+  )  by  CO16  produces  CO(a  3ri)  with  no  re¬ 
ported  observations  of  a  CO(vJ)  product  channel. 

The  energy  transfer  process  reported  here  represents  a 
previously  unobserved  quenching  reaction  of  CO  with  meta¬ 
stable  nitrogen.  We  will  show  that  consideration  of  kinetic 
and  energetic  constraints  identifies  N2(a'  '2~ )  as  the  re¬ 
sponsible  agent.  We  will  also  present  and  discuss  the  results 
of  surprisal  analysis  of  the  observed  vibrational  distributions 
which  indicates  the  rotationally  excited  component  to  be 
formed  by  two  distinct  mechanisms. 

EXPERIMENTS 

These  experiments  were  performed  in  the  COCHISE 
(COld  CHemiexcitation  Infrared  Stimulation  Experiment) 
cryogenic  discharge  afterglow  apparatus  which  is  described 
in  detail  elsewhere. 17  Excitation  of  nitrogen  is  achieved  with 
four  parallel  microwave  discharges  (2450  MHz,  50  W)  of 
flowing  N2/Ar  mixtures  at  ~  1  Torr  total  pressure.  A  dia¬ 
gram  of  the  reaction  chamber  is  shown  in  Fig.  1.  After  exit¬ 
ing  the  discharge  tubes,  the  gas  expands  into  a  low  pressure 
(~3  mTorr),  cryogenicaily  pumped  chamber  (—20  K), 
where  the  molecules  enter  the  collimated  field  of  view  of  a 
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FlO.  I.  Digram  of  tlie  COCHISE  reaction  dumber.  The  phyncal  dunen- 
aoasof  thecdl  are  0.6  at  in  length  and  0.4  ai  m  diameter. 


«««nmg  monochromator/infrared  detector  assembly.  Resi¬ 
dence  times  in  the  discharge  tubes  are  on  the  order  of  3  to  S 
ms;  an  average  time  of  flight  of  0.5  ±0.1  ms  is  required  for 
the  gases  to  exit  the  discharge  tubes  and  enter  the  field  of 
view.  Opposing  flows  of  argon/carbon  monoxide  are  used  to 
create  a  quasistatic  interaction  region  along  the  centerline  of 
the  field  of  view,  resulting  in  partial  rethermalization  of  the 
expansion  coded  rotational  distributions.  The  gas  residence 
time  in  the  field  of  view  is  ~0.3  ms.  Gaseous  helium  refriger¬ 
ant  maintain*  all  internal  temperatures  at  20  K,  excepting 
the  gas  lines  and  optics  which  are  maintained  by  resistive 
heating  dements  at  80  and  40  K,  respectively.  The  tempera¬ 
ture  of  the  bath  gas  in  the  reaction  cell  is  therefore  80  K. 

The  infrared  emissions  are  observed  by  a  cryogenic  0.3 
m  Czemy-Tumer  monochromator  equipped  with  a  liquid- 
helium-cooled  arsenic-doped  silicon  detector  and  a  grating 
Mazed  at  3/rm.  A  chopper  located  in  front  of  the  monochro¬ 
mator  entrance  slit  modulates  the  signal  at  23  Hz.  Data  col¬ 
lection  is  performed  with  a  computer-interfaced  lock-in  am¬ 
plifier.  The  absolute  uncertainty  in  the  wavelengths  (due  to 
monochromator  drive  error )  of  the  data  is  ±  0.003 /im.  The 
data  were  corrected  for  instrument  responsivity  using  black- 
body  calibration  spectra  taken  in  the  300-370  K  range.  The 
absolute  uncertainty  in  the  accuracy  of  the  Mackbody  tem¬ 
perature  is  ±  3  K,  which  results  in  a  relative  error  of 
±  14%  for  4.0/6. 0^m  intensity  ratios.  Spectra  were  taken 
for  Nj/Ar  mixtures  with  N2  mole  fractions  of 0.005  to  0. 12, 
with  a  mass-balanced  counterflow  of  CO/A r  with  CO  mole 
fractions  of  0.018  to  0.33.  The  data  were  taken  at  a  spectral 
resolution  typically  of  0.013  ftm  [full  width  at  half-maxi¬ 
mum  (FWHM)J. 

RESULTS 

In  all,  24  emission  spectra  of  the  CO  fundamental  region 
were  obtained  at  various  N2  and  CO  mole  fractions.  The 
conditions  are  timilar  to  those  in  which  nitric  oxide  chemilu¬ 
minescence  was  examined.1*  The  general  features  of  the 
emissions  are  relatively  invariant  with  CO  mole  fraction,  but 
the  intensity  of  the  CO(  1-0)  emission  exhibits  a  strong  N2 
mole  fraction  dependence  as  shown  in  Fig.  2.  The  relative 


PIG.  2.  Data  comparison.  The  solid  line  shows  a  typical  spectrum  obtained 
under  condition  of  12%  discharged  N,/Ar  reacting  with  a  coun  inflow  of 
33%  CO/ At.  The  dashed  line  shows  data  taken  at  3.17%  N,/Ar  with  iden¬ 
tical  counterflow  conditions.  Shown  below  are  the  band  centers  for  all  of  the 
observed  transitions. 


intensity  of  the  CO(  1-0)  band  changes  by  a  factor  of  4  for  a 
factor  of  4  change  in  the  N2  mole  fraction. 

The  envelope  degrading  to  the  red  of  the  CO(  1-0)  band 
center  at  4.666  /an  is  due  to  the  Av  =  1  progression  from 
v<  14.  The  five  sharp  red-degraded  features  to  the  blue  of  the 
CO(  1-0)  band  center,  which  have  a  spacing  of  29.2  ±  1.4 
cm”  ’,  cannot  be  attributed  to  CO  vibrational  emission  with 
a  rotational  distribution  the  same  as  the  bath  gas  (80  K). 
These  bands  have  not  been  previously  observed  in  published 
spectra  of  CO  fundamental  emission  taken  at  higher  pres¬ 
sures  ( >  1  Torr). 19-21  At  low  pressure  with  the  COCHISE 
apparatus,  these  features  are  observed  under  all  conditions 
that  produce  the  CO  fundamental  emission.  They  exhibit  no 
additional  structure  even  at  the  highest  resolution  employed 
(0.0067 /um).  Undo*  low  nitrogen  mole  fraction  conditions, 
in  which  the  CO(  1-0)  emission  intensity  is  greatly  reduced, 
eight  bands  are  observed  with  three  progressing  into  the 
CO(  Au  =  1)  envelope.  This  is  a  lower  limit  considering  the 
overlap  of  the  features  to  the  red  of  4.6 /im. 

These  emissions  are  not  observed  in  the  absence  of  nitro¬ 
gen  in  the  discharge  mixture.  Thus,  argon  metastables  and 
residual  ions  do  not  contribute  to  the  excitation.  These  fea¬ 
tures  do  not  correspond  to  CN  ( v),  electronic  N2  transitions, 
or  NCO.  The  spectral  shape  of  these  features,  sharply 
peaked  and  degraded  to  the  red,  is  similar  to  the  nitric  oxide 
A-branch  bandheads  identified  in  spectra  of  chemilumines¬ 
cence  produced  from  the  reaction  of  discharged  nitrogen 
with  oxygen.  '*  CO  forms  bandheads  similarly;  at  sufficiently 
high  rotational  excitation,  bandheads  will  form  in  the  R 
branches  and  the  P  branches  will  extend  to  the  red  and  do 
not  form  bandheads.  We  have  used  a  spectral  generation 
technique  to  predict  the  band  shapes  and  positions  of  the  CO 
/(-branch  bandheads.  The  methodology  and  results  will  be 
presented  in  detail  in  the  following  section.  Using  this  tech¬ 
nique  the  sharp  red-degraded  features  have  been  positively 
identified  as  CO  /(-branch  bandheads.  These  features  are 
sufficiently  intense  that  we  have  been  able  to  determine  an 
optimum  rotational  distribution.  Additionally,  the  P 
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branches  of  the  rotationally  excited  component  associated 
with  the  A-branch  band  heads  are  identifiable.  This  verifies 
the  spectral  assignment  and  confirms  the  appropriate  choice 
of  rotational  distribution. 

SPECTRAL  ANALYSIS 

The  data  have  been  analyzed  using  a  spectral  generation 
linear  least-squares  fitting  technique.22  A  computed  infinite 
resolution  spectrum  is  convolved  with  the  instrument  scan 
function  ( in  this  case,  a  symmetric  triangle  with  full  width  at 
half-maximum  as  the  spectral  resolution )  to  create  simulat¬ 
ed  spectra  for  each  vibrational  transition.  The  simulated 
spectra  are  then  fit  to  each  experimental  spectrum  using  a 
linear  least-squares  method  yielding  a  determination  of  the 
product  of  the  upper  state  density  and  the  spontaneous  emis¬ 
sion  coefficient  of  the  transition  N„  A^  P- . 

The  spectroscopic  data  used  in  this  study  are  from 
Huber  and  Herzberg.21  These  data  were  sufficient  to  repro¬ 
duce  the  line  positions  of  both  the  rotationally  relaxed  emis¬ 
sion  features  and  the  A-branch  bandheads  adequately. 

The  rotational  temperature  used  to  reproduce  the  CO 
vibrational  progression  to  the  red  of  4.6  fim  was  determined 
empirically  by  reproducing  the  observed  branch  structure. 
The  optimum  temperature  was  determined  to  be  80  K, 
which  is  the  same  temperature  as  the  bath  gas.  Thus,  this 
emission  system  is  rotationally  thermalized.  The  rotational 
distribution  was  treated  by  a  simple  Boltzmann  expression, 
so  that  band-integrated  vibrational  number  densities  and 
transition  probabilities  were  used.  The  transition  probabili¬ 
ties  were  calculated  using  the  dipole  moment  function  of 
Chackerian  et  at.24  Fourteen  vibrational  levels  of  80  K  CO 
emission  have  been  positively  identified  from  the  data.  This 
corresponds  to  vibrational  excitation  of  3.42S  eV. 

The  band-averaged  transition  probabilities  used  for  the 
rotationally  excited  component  are  the  same  as  those  for  the 
80  K  CO  emission.  We  believe  these  values  to  be  accurate 
since  the  CO  ground  state  is  '£  and  therefore  not  subject  to 
spin  uncoupling  at  high  J'  as  discussed  for  NO.1*  Thus, 
Honl-London  scaling  of  the  band-averaged  transition  pro¬ 
babilities  should  be  sufficient  to  determine  the  populations 
of  the  rotationally  excited  component  accurately. 

The  0.013  fim  resolution  (3  cm-1  at  3.0  fim)  of  the 
spectral  data  is  insufficient  for  rotational  resolution  of  the  A- 
branch  bandheads  so  the  chosen  rotational  distributions  are 
those  which  best  reproduce  the  band  shapes  and  peak  posi¬ 
tions.  To  fit  the  A-branch  bandhead  features  adequately  us¬ 
ing  a  Boltzmann  distribution  requires  temperatures  of  ap¬ 
proximately  20000  K.  Although  the  fits  to  the  lower 
vibrational  levels  are  adequate,  higher  vibrational  levels  are 
poorly  fit  and  are  better  described  by  lower  Boltzmann  rota¬ 
tional  temperatures.  This  is  evidence  for  an  anticorrelation 
between  rotational  and  vibrational  excitations;  i.e.,  the  low¬ 
est  vibrational  levels  contain  the  most  rotational  excitation. 
Such  anticorrelations  have  been  observed  in  photolysis  ex¬ 
periments,  abstraction  reactions,  and  energy  transfer  pro¬ 
cesses.1'*4  In  these  instances,  the  observed  rotational  distri¬ 
butions  are  sometimes  best  described  by  a  statistical  model.3 
Such  a  model  distributes  the  population  statistically  over  all 
accessible  states.  The  model  that  incorporates  the  observed 


anticorrelation  between  vibrational  and  rotational  excita¬ 
tion  is  gi  ven  by 

A»(y)«(2y-H){(£r-A„)-£,},/2.  (3) 

where  ET  is  the  total  energy  available  for  product  states,  E, 
is  the  vibrational  energy,  and  Ej  is  the  energy  of  the  particu¬ 
lar  rotational  level.  Figure  3  contrasts  the  relative  popula¬ 
tion  distributions  predicted  by  Boltzmann  and  statistical 
models.  The  statistical  model  contains  greater  relative  popu¬ 
lation  in  the  higher  rotational  levels  at  the  expense  of  the 
lower.  Thus,  rotational  bandhead  formation,  which  requires 
significant  population  of  rotational  levels  in  the  region  of  the 
Fortrat  parabola  vertex  =  92  for  v  =  1 ),  is  readily 
facilitated  by  a  statistical  distribution. 

We  have  examined  several  values  of  Er  to  determine 
which  best  reproduces  the  positions  and  spectral  shapes  of 
the  observed  eight  bandheads.  Values  of  E1  below  3.0  eV 
reproduce  the  lower  vibrational  levels  well,  but  the  predicted 
bandheads  at  higher  vibrational  levels  are  broadened  and  red 
shifted.  The  responsible  mechanism  may  be  seen  from  Fig.  3. 
At  higher  £„  and  lower  Ej,  the  relative  population  of  the 
rotational  levels  near  the  region  of  the  reversal  decreases 
which  causes  the  observed  effects.  Values  of  Er  greater  than 
4.0  eV  predict  bandhead  formation  for  vibrational  levels  of 
ten  or  greater,  which  have  not  been  unambiguously  identi¬ 
fied  in  the  data.  We  prefer  a  value  of  Er  between  these  two 
bounds.  Comparing  fits  using  several  values  within  this 
range,  we  have  determined  Er  =  3.7  eV  to  be  optimum.  Val¬ 
ues  within  ±  0.2  eV  of  the  optimum  value  produce  adequate 
fits  with  only  small  differences. 


J-  80-100 


FIG.  3.  Comparison  of  various  rotational  distributions.  Shown  are  the  rela¬ 
tive  rotational  distributions  for  a  10*  K  Boltzmann  and  several  statistical 
distributions. 
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TAMlBwplta  of  CO  h— rtfcwrt  fcr—tton  E,  -  3.70  eV. 
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Table  I  shows  the  energetics  of  CO  bandhead  formation 
for  Er  of  3.7  eV.  Shown  are  the  values  for  EK,  (ET  —  E„), 
the  rotational  level  corresponding  to  the  vertex  of  the  For- 
trat  parabola,  the  rotational  energy  corresponding  to  this 
value  of/,  and  the  maximum  rotational  level  allowable  from 

E* 

Figure  4  shows  a  typical  fit  to  the  data  using  14  vibra¬ 
tional  levels  of  80  K  CO  emission  and  eight  vibrational  levels 
with  a  statistical  rotational  distribution  corresponding  to  Er 
—  3.7  eV.  All  of  the  principal  spectral  features  are  accurate¬ 
ly  reproduced.  The  only  features  not  entirely  reproduced  fall 
within  the  5. 7-6.5  fim  region.  Figure  5  shows  an  enlarged 
view  of  this  region  from  Fig.  4.  The  fit  is  shown  by  the  heavy 
line.  The  spacing  of  the  bands  shown  in  Fig.  5  is  ~6  cm" '. 
These  features  have  been  observed  in  all  spectra  containing 
sufficient  intensity  in  this  wavelength  region.  The  fit  shows 
excellent  reproduction  of  the  spectral  shapes,  but  does  not 
match  the  absolute  intensity.  Figure  4  shows  that  the  80  K 
CO  emission  does  not  contribute  to  this  spectral  region;  the 
discrete  features  are  reproduced  by  the  P  branches  of  the 
rotationally  excited  component.  Inclusion  of  v  >  8  of  the  ro- 
tationally  excited  component  improves  the  fits  in  this  wave¬ 
length  region.  Due  to  overlap  with  the  CO  envelope,  how¬ 
ever,  unique  determination  of  the  populations  is  not  possible, 
so  we  have  generally  excluded  these  bands  from  the  fits.  By 
comparing  fits  with  different  vibrational  contributions  to  the 
rotationally  excited  component,  we  have  determined  that 
the  majority  of  the  intensity  in  the  5.7-6. 5  fim  region  arises 
from  the  higher  vibrational  levels  ( v '  =  6-8).  The  relative 
line  spacings  may  be  calculated  (omitting  D,,  correction) 
from 

Av  =  />(/)-/>(/  +  1 )  =  (2?„  +  B„. ) 

-  (ff„  -BV-H2J+  1). 

(4) 

For  v'  =  6-8  this  wavelength  region  contains  F-branch  line 
spacings  of  —6  cm-1  for  rotational  levels  50-80.  The 
N2(  v'  =  1  —  B  ’nf,  v"  =  0)  emission  in  the  6.3-6.6 


4.0  4.S  5.0  5.5  6.0  6.5 


WAVELENGTH  (pm) 

FIG.  4.  Data  (light  line)  and  best  fit  (dark  line)  (a)  to  80  K  CO  fundamen¬ 
tal  emission  (#'  =  1-14)  and  CO  rotationally  elicited  bands  using  a  statisti¬ 
cal  distribution  with  E,  =  3.7  eV.  The  data  are  the  same  as  used  for  the 
dashed  line  in  Fig.  2.  The  spectral  resolution  is  0.01 3//m.  The  80  K  CO  and 
rotationally  excited  simulated  spectra  which  comprise  the  best  fit  in  (a)  are 
shown  in  (b)  and  (c),  respectively. 


fi m  region,  which  wc  have  previously  examined  and  report¬ 
ed,22  and  signal-to-noise  considerations  impede  determina¬ 
tion  of  the  full  extent  of  the  F-branch  structure.  The  repro¬ 
duction  of  these  features,  however,  by  fits  to  the  rotationally 
excited  component  confirms  the  identification  of  the  CO  R- 
branch  bandheads. 


FIG.  5.  Enlarged  view  of  the  5.7  to  6.5  fim  region  of  Fig.  4.  The  data  is 
shown  by  a  light  line,  the  fit  by  a  dark  line.  The  data  shown  here  exhibit  a 
small  wavelength  offset  from  the  fit  well  within  the  monochromator  drive 
error. 
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AH  24  spectra  have  been  fit  using  14  80  K  CO  emission 
bands  and  eight  vibrational  leveh  of  the  rotationally  excited 
component,  as  shown  in  Fig.  4.  The  statistical  rotations!  dis¬ 
tributions  employed  in  these  fits  provide  a  better  reproduc¬ 
tion  of  the  data,  both  in  the  A-branch  spectral  band  shapes 
and  the  absolute  intensity  of  the  P  branches,  than  do  Boltz¬ 
mann  rotational  distributions.  We  consider  the  statistical 
model  employed  here  to  be  a  more  accurate  representation  of 
the  true  rotational  distributions,  but  it  cannot  be  considered 
to  be  a  unique  determination.  Signal-to-noise  considerations 
and  band  overlap  do  not  permit  unambiguous  identification 
of  the  vibrational-level  dependent  rotational  distributions. 
The  true  distributions  may  deviate  from  a  purely  statistical 
model,  possibly  containing  some  Gaussian  character.  How¬ 
ever,  the  model  employed  here  has  successfully  demonstrat¬ 
ed  an  anticorrelation  between  vibrational  and  rotational  ex¬ 
citation  and  permitted  a  bound  of  ~  3.7  eV  to  be  estimated 
for  the  rotationally  excited  component.  Determination  of 
this  bound  and  its  similarity  to  the  maximum  vibrational 
excitation  of  the  80  K  CO  component  3.42$  eV  suggests  that 
these  two  components  arise  from  the  same  excitation  pro¬ 
cess. 

KINETIC  INTERPRETATIONS 

The  kinetics  of  processes  occurring  in  the  COCHISE 
reaction  chamber  have  been  described  previously.18  Owing 
to  the  low  number  densities  in  the  reaction  zone,  the  short 
residence  time  in  the  field  of  view  (0.3  ms)  and  the  long 
radiative  lifetimes  for  the  observed  infrared  chemilumines¬ 
cent  processes,  reactions  and  quenching  (excluding  rota¬ 
tional  )  of  the  excited  species  created  in  the  reaction  zone  can 
be  neglected.  Thus,  vibrational  quenching  and  CO(o)  up 
pumping  [reaction  (lb)]  can  be  ignored.  The  kinetics  of 
[CO(t;,/)  ]  are  therefore  in  the  steady  state  according  to  the 
production  rate  and  the  lifetime  of  the  excited  species  in  the 
field  of  view 

d  [CO (vj)\/dt  =  k  [Af  •]  [CO]  -  r- 1 '(CO(M)  ]  =  0, 

(5) 

where  Af*  denotes  the  excited  species  created  in  the  micro- 
wave  discharges  responsible  for  the  observed  excitation,  k  is 
the  excitation  rate  coefficient,  and  rm  is  the  residence  time 
in  the  field  of  view,  0.3  ms. 

Figure  6  shows  the  absolute  vibrational  populations  for 
the  fit  in  Fig.  4.  The  populations  for  the  rotationally  excited 
component  exceed  those  of  the  80  K  component  at  all  vibra¬ 
tional  levels  except  v=  1,2. 

As  illustrated  by  Fig.  2  the  lowest  80  K  CO  vibrational 
populations  exhibit  a  dependence  on  N2  mole  fraction.  We 
have  determined  previously"1  that  only  N2(i>)  exhibits  a 
strong  variation  in  its  discharge  production  rate  with  N2 
mole  fraction  over  the  range  used  here.  Thus,  relative  in¬ 
creases  in  low  CO  vibrational  population  at  higher  nitrogen 
mole  fraction  must  be  due  to  the  near-resonant  energy  trans¬ 
fer  from  N2(t>)  [reaction  ( la)  ].  Examination  of  the  popula¬ 
tions  for  all  the  spectra  indicates  that  only  CO(v  =  1,2)  are 
affected  by  nitrogen  mole  fraction.  Since  CO(t>)  up  pumping 
cannot  occur  under  our  experimental  conditions,  the  excita¬ 
tion  process  must  be  single-  and  two-quantum  transfer 


FIG.  6.  Determined  population  distribution  vs  E,„  for  the  fit  shown  in  Fig. 
4.  The  thermalized  component  is  shown  as  (O),  with  contribution  to  ther- 
malized  CO(i>  =  1,2)  from  N,(v)  energy  transfer  shown  as  (•).  and  the 
rotationally  excited  component  is  represented  as  (□). 


N2(u)  +  CO—N2(o  —  1)  +  CO(i>  =  1),  (6a' 

N2(p)  +  CO— N2(i>  —  2)  +  CO(u  =  2).  (6b) 

To  determine  the  relative  contribution  of  the  N2(p)  energy 
transfer  process  to  the  CO(v  =  1,2)  populations,  we  have 
plotted  the  ratios  of  these  populations  to  [CO(o  =  3)  ]  as  a 
function  of  nitrogen  mole  fraction.  Figure  7  shows  one  of 
these  plots.  Extrapolation  of  these  curves  to  zero  nitrogen 
mole  fraction  determines  the  relative  CO(i>  =1,2)  popula¬ 
tions  which  arise  from  sources  other  than  energy  exchange 
from  N2(u  =  1,2).  For  the  data  in  Fig.  7,  the  multiplicative 
factors  are  3.0  and  1.5  for  CO (u  —  1,2),  respectively.  The 
multiplicative  factors  for  all  CO  mole  fractions  examined 
have  been  determined  to  be  approximately  the  same. 

Using  these  factors,  the  contributions  to  the  80  K 
CO(v=  1,2)  populations  from  N2(u)  transfer  and  E-V 
transfer  may  be  separated.  The  dotted  line  in  Fig.  6  shows 
the  CO(u=  1,2)  number  densities  corrected  for  N2(u) 
transfer.  The  vibrational  populations  of  the  rotationally  ex¬ 
cited  component  exceed  those  of  the  80  K  component  by 
typically  —30%.  This  ratio  is  similar  for  all  the  analyzed 
spectra,  showing  thereby  that  neither  component  exhibits  a 
dependence  on  nitrogen  or  carbon  monoxide  mole  fraction. 

The  CO(u  =  2)  population  attributable  to  two  quan¬ 
tum  transfer  from  N2(i>)  is  only  6%  ±  1%  of  that  due  to 
single  quantum  transfer.  This  ratio  is  related  to  the  rate  coef¬ 
ficients  for  reactions  ( 6a )  and  ( 6b )  in  the  following  manner: 
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DISCHARGE  N2  MOLE  FRACTION 

FIG.  7.  The  plot  of  the  ratio  of  the  [CO(v  =  1,2)1  to  |CO(v=  3)]  popula¬ 
tions  as  a  function  of  nitrogen  mole  fraction.  The  CO  mole  fraction  for  these 
data  was  0.33. 


[CO(tt  =  2)1  £A(Ap  =  2)[N2(i;)] 

[CO(p  =  1)1  2„/t„(Au=  1)[N2(p)1  ' 

Equation  (7)  may  be  resolved  if  the  N2(t>)  distribution  can 
be  determined.  Discharge-flow  measurements  by  Piper  and 
Marinelli  using  Penning  ionization  spectroscopy  indicate 
the  effluent  of  microwave  discharges  employing  the  condi¬ 
tions  encountered  here  may  be  represented  by  a  ~6000  K 
“modified  Treanor”  distribution  in  N2(o).2’  A  Treanor  dis¬ 
tribution  contains  enhanced  populations  at  higher  vibration¬ 
al  levels,  relative  to  a  Boltzmann  distribution,  that  are  creat¬ 
ed  by  collisional  up  pumping.14  Using  the  modified  Treanor 
distribution,  we  have  determined  k„(Ay  =  2)/fc„(Ai>  =  1) 
=  0. 10  ±  0.04.  Owing  to  the  known  increase  in  ku  as  a  func¬ 
tion  of  i>,IUI  the  ratio  determined  here  is  likely  representa¬ 
tive  of  higher  v  ( probably  7  to  8 ) .  This  value  falls  well  within 
the  limits  of  0.2  to  0.025  found  for  all  vibrational  levels.26 

Figure  8  shows  the  average  population  distributions  for 
the  80  K  and  rotationally  excited  components,  normalized 
separately  and  shown  with  one  standard  deviation  error 
bars.  The  population  distributions  from  1 1  spectra  were  cho¬ 
sen  for  this  average.  The  spectra  containing  large  contribu¬ 
tions  to  80  K  CO(v  =  1,2)  levels  have  been  excluded  since 
this  emission  envelope  overlaps  the  v  =  5-8  levels  of  the  ro¬ 
tationally  excited  component  interfering  with  reliable  popu¬ 
lation  determination. 

The  distributions  shown  in  Fig.  8  are  similar,  but  the 
distribution  of  the  rotationally  excited  component  appears 
to  be  relatively  flat  above  v  =  4.  We  have  examined  tiie  data 
carefully  and  determined  this  trend  to  be  accurate.  Fits  to 
the  data  using  a  fixed  relative  vibrational  distribution  of  the 
rotationally  excited  component  following  the  distribution  of 
the  80  K  component  seriously  underfits  the  data  in  the  re- 


Ex103cm-i 

FIG.  8.  Avenged  relative  population  distributions  for  1 1  spectra.  Both  80 
K  (•)  CO  and  the  rotationally  excited  components  (□)  are  shown,  but 
have  been  normalized  independently.  The  error  bars  represent  one  standard 
deviation. 


gionoft/  =  5-8  ofthe /{-branch  bandheads  and  in  the  region 
of  the  P  branches. 

EXCITATION  MECHANISM 

One  possible  mechanism  for  rotational  excitation  of  CO 
is  energy  transfer  from  translationally  excited  atoms  ema¬ 
nating  from  the  discharge.  Translationally  hot  H  atoms  have 
been  reported  to  excite  VJt  states  of  CO.27  We  may  discount 
such  processes  in  our  apparatus  since  we  can  find  no  mecha¬ 
nism  for  hot  atom  formation  and  any  “hot”  atoms  that  could 
form  in  the  1  Torr  discharge  tubes  would  be  rapidly  accom¬ 
modated.  The  number  density  of  such  species  reaching  the 
field  of  view  must  therefore  be  miniscule. 

Collisional  quenching  of  high  CO  vibrational  levels  into 
high  rotational  levels  of  lower  vibrational  levels  ( V,R  trans¬ 
fer),  analogous  to  processes  observed  for  HF,28  29  may  also 
be  dismissed.  To  excite  the  high  J '  levels  we  observed  would 
require  nearly  gas  kinetic  multiquantum  (Ai><8)  quenching 
ofCO(u).  This  is  unlikely  since  direct  measurements  of  the 
p-dependent  quenching  of  CO(p)  by  CO, 30  indicate  the  to¬ 
tal  quenching  rate  coefficients  to  be  less  than  1  %  gas  kinetic. 
Under  quasiresonance  conditions,  collisional  quenching  of 
diatomics  in  low  v,  high  J  levels  into  high  p,  low  J  levels  (R,V 
transfer)  may  have  rate  coefficients  of  10”  "  cm’ s-1.  ”  The 
resonance  conditions  are  given  by 

oi„  / (Oj  =  at,  MB,J.  ( 8 ) 

For  CO,  resonance  is  satisfied  at  J '  ~  1 20  which  is  populated 
only  at  the  lowest  vibrational  level  for  the  statistical  model. 
Therefore,  R,V  transfer  should  not  contribute  significantly 
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to  the  higher  vibrational  levels  of  the  observed  thermalized 

With  other  possibilities  excluded,  the  mechanism  for 
CO(uJ)  excitation  must  be  energy  transfer  from  one  or 
more  of  the  metastable  nitrogen  species  created  in  the  dis¬ 
charge.  We  have  previously  determined  that  to  account  for 
the  CO(vJ)  excitation  observed  here,  the  energy  transfer 
reaction  must  be  at  least  ~  3.7  eV  exoergic.  We  use  this  con¬ 
straint,  combined  with  kinetic  considerations  to  identify 
Nj(a*  '2  ~  )  as  the  species  responsible  for  the  observed  CO 
excitation. 

Summing  the  populations  of  the  80  K  and  rotationally 
excited  components,  a  product  k[A/  *]  ~0.5  s- 1  is  required 
to  account  for  the  observed  emissions.  The  measured 
quenching  kinetics12-36  for  many  of  the  metastable  nitrogen 
species  are  shown  in  Table  II.  This  table  contrasts  the  num¬ 
ber  densities  of  the  metastable  species  required  in  the  inter¬ 
action  zone  to  account  for  the  observed  emissions  with  the 
number  densities  determined  from  modeling  studies  or  di¬ 
rect  determination.  '*  The  number  density  calculations  have 
employed  the  room  temperature  rate  coefficients  in  the  ab¬ 
sence  of  data  at  80  K. 

Quenching  reactions  of  metastable  nitrogen  atoms  are 
shown  in  reactions  (9)  and  ( 10): 

N(2Z>)  +  CO— N(4S)  +  CO(o<9),  A£  =  2.38  eV,  (9) 
N(2/,)+CO-.N(2Z>)+CO(u<4),  AE  =  1.196  eV, 

(10a) 

-N(4S)+CO(u<15),  A£  «  3.576  eV 

(10b) 

N(2Z>)  may  deposit  2.38  eV  into  CO  ro vibrational  states, 
well  below  the  —  3.7  eV  needed  to  account  for  the  observed 
emissions.  Additionally,  the  kinetics  of  this  reaction  indicate 
that  contributions  to  the  spectra  from  this  source  would  be 
negligible.  The  energetics  of  N(2/*)  quenching  to  form 
N(4S)  matches  the  required  ~3.7  eV  well.  However,  the 
recently  determined  rate  coefficient32  for  quenching  of 
N(2P)  by  CO  is  several  orders  of  magnitude  too  small  for 
reactions  (10a)  and  (10b)  to  be  the  CO(M)  excitation 
mechanism. 


TABLE  II.  Measured  CO  quenching  kinetics  of  NJ  JV  •. 


Number  density  in  COCHISE  (cm'-’) 


Species 

K0  (300  K) 
(cm"’  s_ ') 

Required  for 
CO(vJ)  excitation 

Estimated  for 
interaction  zone 

N  CP) 

<l.5< 

-  14)' 

>303) 

1-3(9) 

N  (JZ» 

l.7( 

-  12)" 

>3(11) 

3-10(9) 

NjM  ’2„‘ ). 
v  =  0 

l.5( 

-  12)b' 

>3(11) 

1-3(9) 

o  =  4 

!.9( 

-  1 1 )** 

3(10) 

1-3(9) 
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2.8< 

-  I0)*-r 

2(9) 

2(6) 

N2(e-  'E.  ) 

IK 

-  ior* 

4(9) 

3(9) 

*  Includes  excitation  of  COM  'll).  'Reference  33. 

*  Includes  excitation  of  CO(o  ’ll ).  •  Reference  36. 

'  Reference  32.  Notation:  1  (  —  10)  =  1 X 10  “ 

"Reference  33. 

'Reference  34. 


Excitation  of  CO(v<14)  by  energy  transfer  from 
N2(.T,i/)  requires  multiquantum  transfer  from  t/<14.  The 
contributions  to  the  80  K  CO(v  =  1,2)  populations  deter¬ 
mined  from  N2(y)  energy  transfer  are  consistent  with  a  pre¬ 
ferred  channel  for  single-quantum  exchange  and  a  less  fa¬ 
vored  channel  for  two-quantum  exchange.  Therefore,  we  do 
not  believe  multiquantum  exchange  from  N2(X,v')  can  be 
the  responsible  excitation  mechanism. 

N2M  32  ~  )  quenches  with  CO  to  form  CO(a  3I1).34 
Although  the 

N}(^32-)  +  co-.N2m  +  co(fl3n),  an 

CO(a3m  +  CO-2CO(Jr,5<u<13)  (12) 

rate  coefficient  for  the  quenching  process  is  known,34  the 
absolute  yield  for  reaction  (11)  has  not  been  measured.  Re¬ 
action  of  CO(a)  with  another  CO  molecule  produces  vibra- 
tionally  excited  CO37-38  with  a  nearly  gas-kinetic  rate  coeffi¬ 
cient  for  the  process  (k,2~-lxl0_l°  cm3  s~ 1  at  room 
temperature). 3,40  Quenching  of  N2M)  may  produce  rovi- 
brationally  excited  CO  directly  from  a  branch  in  reaction 
(11).  Since  N2M)  lies  6.17  eV  above  N2(A>  =  0),  such  a 
process  would  be  sufficiently  energetic.  The  N2M )  quench¬ 
ing  rate  coefficient,  however,  is  too  small  for  this  process  or 
any  derived  reactions,  such  as  reaction  ( 12),  to  be  responsi¬ 
ble.  Additionally,  the  activation  barrier  determined  for  this 
process  by  Slanger  et  a  l 41  indicates  the  rate  constant  may  be 
smaller,  <  10“ 13  cm3  s~ ',  at  80  K. 

The  N2(  W3bu,w  1  A„ )  states  he  7.36  and  7.35  eV  above 
the  ground  state,  respectively,  and  are  sufficiently  energetic 
to  produce  the  observed  CO(  vj)  emissions.  The  quenching 
reactions  of  these  two  species  have  not  been  reported 
in  the  literature,  however.  We  have  determined 
N2(IF3Ab,i;=  1-5)  andN2(u>  'A„,u  =  0-2)  number  densi¬ 
ties  in  COCHISE  directly  from  their  IR  radiance  over  the  2- 
4/xm  region.22  Their  concentrations  in  the  interaction  region 
are  typically  3x10*  and  1 X 10*  molecules  cm  " 3,  respective¬ 
ly.  Even  if  these  species  quench  CO  with  rate  coefficients 
near  gas  kinetic,  they  cannot  account  for  the  observed 
CO(u^f)  emissions. 

Quenching  of  N2(a  'll)  is  rapid,  forming 
COM  'n).15'35  Although  the 

n2(o  'id  +  co-N2(;r,io  +  com  'n)  03) 

yield  of  COM)  formation  from  this  reaction  is  large,15  a 
channel  forming  CO  (vj)  is  possible.  N2(a)  may  deposit 
any  fraction  of  the  available  8.4  eV  into  CO  rovibrational 
states  from  such  a  process.  However,  the  radiative  lifetime  of 
N 2(a)  is  short,  56  ±  4/is,42  which  makes  the  number  den¬ 
sity  of  this  species  in  the  interaction  zone  to  be  too  small  to 
account  for  the  observed  emissions.  The  upper  limit  for 
N2(o)  number  density  shown  in  Table  II  has  been  deter¬ 
mined  from  the  noise  level  of  discharged  Ar/N2  spectra  at 
the  wavelength  N2(o  —  a')  features  would  occur,22  and  is 
consistent  with  the  expected  radiative  decay  of  N2(j)  in  the 
0.5  ms  collisionless  expansion  between  the  discharge  exit 
and  the  reaction  zone. 

The  possible  quenching  reactions  of  the  N2(a’  '2  ~  ) 
with  CO  are 
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Nj (a'  '2~  )  +  CO-»N2(2’,ir)  +  C 0(A  'n),  ( 14a) 

-N2M,t>)  +  CO(X,vJ),  ( 14b) 

-NjM.u)  +CO(/'2-,2)'A). 

(14c) 

Nj(a')  may  deposit  up  to  8.5  eV  into  CO  electronic,  vibra¬ 
tional,  and  rotational  states.  The  branching  ratio  forming 
COM)  nas  been  measured  as  20  This  value,  how¬ 

ever,  was  based  on  an  80  fts  lifetime  for  the  N 2(a)  state 
which  has  recently  been  revised  to  be  56  ±  4  ft s.4J  The  cor¬ 
rected  COf  A )  branching  ratio  is  30%  ±8%.  The  remaining 
fraction  must  be  divided  between  the  other  energetically 
accessible  spin-allowed  channels.  These  channels  are  rovi- 
brationally  excited  COM)  [reaction  ( 14b)  ]  and  I '2~  and 
D  'A  states  [  reaction  ( 14c)  ] .  At  room  temperature,  the  to¬ 
tal  quenching  rate  constant  for  reaction  ( 14)  is  1. 1 X  10'  1,1 
cm3  s~  Table  II  shows  that  the  kinetics  of  this  reaction 
are  sufficient  to  account  for  the  observed  emissions. 

Since  N2(a')  is  the  only  species  present  in  the  interac¬ 
tion  zone  that  satisfies  both  the  kinetic  and  energetic  con¬ 
straints,  we  postulate  that  reaction  ( 14b)  must  be  responsi¬ 
ble  for  the  observed  CO(vJ)  excitation.  For  this  hypothesis 
to  be  correct,  the  branching  fraction  for  reaction  ( 14b )  must 
be  relatively  large,  constraining  reaction  (14c)  to  be  small 
and,  owing  to  the  short  radiative  lifetime  of  COM )  •  the  80  K 
populations  must  reflect  a  contribution  from  COM)  radia¬ 
tive  cascade. 

DISCUSSION 

In  the  preceding  discussion,  we  identified  N2(a'  ’1~  ) 
to  be  the  species  responsible  for  the  observed  CO  (tv/)  exci¬ 
tation.  We  show  in  this  section  that  simple  models  incorpor¬ 
ating  reactions  ( 14a)  and  ( 14b)  can  account  for  the  ob¬ 
served  vibrational  distributions  for  both  the  80  K  and 
rotationally  excited  components.  Surprisal  analysis  of  the 
rotationally  excited  component  vibrational  distribution  in¬ 
dicates  two  mechanisms  are  responsible  and  vibrational  lev¬ 
els  greater  than  eight  are  predicted.  We  then  demonstrate 
that  the  80  K  vibrational  distributions  can  be  modeled  by 
contributions  from  COM)  radiative  cascade  and  the  rota¬ 
tionally  excited  component.  The  latter  may  arise  either  from 
rotational  relaxation  or  initial  bimodality  in  the  energy 
transfer  reaction. 

Surprisal  theory43**  postulates  the  existence  of  an  expo¬ 
nential  gap  law  for  an  individual  metathetic  reaction 

P„  =  P°(v)exp(  -  Af„)/exp(A„),  (15) 

where/,  is  the  fraction  of  the  reaction  exoergicity  appearing 
as  vibrational  energy  in  the  product,  P(v)  is  the  observed 
relative  vibrational  population,  and  P°(v )  is  the  statistical 
“prior”  distribution  obtained  when  all  final  translational, 
rotational,  and  vibrational  states  are  equally  probable.  Thus, 
the  ratios  of  the  observed  and  statistical  vibrational  popula¬ 
tions  should  be  exponential  in  the  vibrational  energy  E,„ 
with  the  exponential  fall-olf  constant  A  quantifying  the  de¬ 
parture  of  the  observed  distribution  from  a  completely  statis¬ 
tical  product  distribution.  While  this  theory  does  not  neces¬ 
sarily  hold  for  all  chemical  interactions,  it  has  proved  useful 
in  the  analysis  of  rovibrational  product  distributions  from 


several  photochemical  and  reactive  interactions. 

Both  the  vibrational  extent  of  the  80  K  CO  emission 
(i><14,  E~  3.4  eV)  and  the  rovibrational  excitation  of  the 
rotationally  excited  component  (ET~  3.7eV)  can  be  attrib¬ 
uted  to  an  energy  transfer  process  of  3. 5-3.7  eV.  Using  a 
prior  distribution  for  3.7  eV  determined  from  the  usual  rela¬ 
tionship  for  a  vibrating  rotator,43-44 

P»  =  (l  -f,)in  £(l-/,)3/J  (16) 

v= 0 

the  vibrational  surprisal  plot  for  the  vibrational  populations 
of  the  rotationally  excited  component  has  been  plotted  in 
Fig.  9.  The  change  in  slope  indicates  two  dynamic  mecha¬ 
nisms  are  responsible  for  the  rotationally  excited  compo¬ 
nent:  one  accounting  principally  for  v  =  1-4  and  another  for 
the  higher  vibrational  levels.  The  relative  populations  of  the 
rotationally  excited  component  may  be  reproduced  by  a 
model  composed  of  two  such  mechanisms  as  shown  in  Fig. 
10.  The  model  is  comprised  by  roughly  equal  contributions 
from  a  low  v  excitation  process  and  a  “statistical”  process 
which  contributes  to  all  vibrational  levels.  The  distribution 
used  for  the  statistical  process  is  the  prior  calculated  for  3.7 
eV.  This  model  predicts  populations  for  v  >  8  of  the  rotation- 
ally  excited  component.  Absolute  populations  of  these  levels 
cannot  be  determined  from  the  spectra  due  to  the  overlap  of 
these  bands  with  the  80  K  CO  emission  features.  However, 
we  have  performed  fits  to  the  data  using  v  =  1-14  of  the 
rotationally  excited  component  following  the  modeled  dis¬ 
tribution  shown  in  Fig.  10.  The  resulting  fits  do  not  exhibit 
any  marked  differences  from  fits  excluding  v  >  8  of  the  rota¬ 
tionally  excited  component  and  the  match  to  the  absolute 
intensity  in  the  P-branch  region  is  generally  improved.  Un¬ 
der  these  conditions,  the  populations  determined  for  t>>  4  of 
the  80  K  component  are  decreased  by  an  average  of  —40%. 


FIG.  9.  Vibrational  surprisal  plot  for  the  rotationally  excited  vibrational 
populations.  The  populations  are  those  from  Fig.  8. 
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FIG.  10.  Model  fit  to  the  rotationally  excited  vibrational  populations.  The 
model,  shown  by  a  solid  line,  has  been  produced  by  equal  contributions 
from  the  k>wi>(v<4)  (— )  mechanism  and  the  statistical  mechanism  (all  v) 
(•). 


Thus,  the  populations  for  t>>8  of  the  80  K  component 
shown  in  Fig.  8  represent  upper  limits. 

Due  to  the  implications  of  the  model  and  the  improve¬ 
ment  in  the  P-branch  fits,  we  believe  that  v>  8  of  the  rota¬ 
tionally  excited  component  are  present  in  the  data.  How¬ 
ever,  higher  resolution  studies  will  be  required  to  confirm 
their  presence  and  determine  the  true  population  distribu¬ 
tion. 

The  80  K  component  must  include  contribution  from 
reaction  ( 14a)  to  the  vibrational  populations  due  to  COM) 
radiative  cascade.  Quenching  of  N2(a'  2.  .v  =  0)  via  reac¬ 
tion  ( 14a)  yields  COM  ln,i><2)/v*  Owing  to  the  short  ra¬ 
diative  lifetime  of  the  COM )  state,  — 10  ns,  all  of  the  COM ) 
created  by  this  process  under  our  experimental  conditions 
will  cascade  radiatively  to  form  CO(tK,9)  with  a  vibrational 
distribution  reflecting  the  known  branching  ratios.  The  rota¬ 
tional  temperature  of  the  COM )  states  formed  in  this  man¬ 
ner  have  been  determined  to  be  approximately  1000  K.*’ 
This  amount  of  rotational  excitation  is  modest  and  the  re¬ 
sulting  COM)  states  formed  by  radiative  cascade  will  be 
rapidly  thermalized  within  (he  residence  time  in  the  field  o. 
view  (20-30  collisions). 

It  is  reasonable  to  expect  that  relaxation  of  the  rotation- 
ally  excited  component  may  also  contribute  to  the  observed 
80  K  component  populations.  Using  only  COM )  radiative 
cascade  and  rotational  relaxation  of  the  rotationally  excited 
component,  we  have  successfully  modeled  the  80  K  vibra¬ 
tional  distribution  as  shown  in  Fig.  1 1 .  The  contribution  to 
the  model  from  the  rotationally  excited  component  is  shown 
in  Fig.  1 1  as  the  low  v  and  statistical  distributions  from  Fig. 
10.  Reproduction  of  the  absolute  distribution  requires 


FIG.  11.  Model  fit  (o  the  80  K  vibrational  populations.  The  model  ( — )  has 
been  produced  by  a  contribution  of  14%  branching  fraction  from  reaction 

( 14a),  ( - )  and  30%  of  the  rotationally  excited  model,  shown  separately 

as  low  v  (--)  and  statistical  contributors  ( — ).  The  relative  populations 
have  been  normalized  to  the  sum  of  the  rotationally  excited  vibrational  pop¬ 
ulations. 


roughly  30%  of  the  rotationally  excited  component  to  be 
thermalized  (relaxed)  and  a  branching  fraction  for  COM) 
excitation  of  14%.  The  determined  contribution  from 
COM )  radiative  cascade  is  about  a  factor  of  2  lower  than 
the  recalculated  branching  ratio  of  Piper.36  Our  result  could 
be  consistent,  however,  if  the  branching  ratio  for  reaction 
(14a)  decreases  somewhat  at  lower  temperatures.  As  an  al¬ 
ternative  to  rotational  relaxation  of  the  rotationally  excited 
component,  the  initial  distribution  may  be  bimodal.  We  can¬ 
not  distinguish,  however,  between  these  scenarios  under  our 
experimental  conditions. 

In  summary,  the  results  of  surprisal  analysis  indicate 
two  dynamic  mechanisms  are  responsible  for  the  rotational- 
ty  excited  component,  one  producing  the  low  vibrational  lev¬ 
els  and  the  other,  a  statistical  mechanism  that  contributes  to 
all  vibrational  levels.  COM )  radiative  cascade  and  relaxa¬ 
tion  of  a  fixed  fraction  of  the  rotationally  excited  component 
are  sufficient  to  account  for  the  80  K  populations. 

The  dynamics  of  the  N,(o')  +  CO  quenching  process 
are  worthy  of  some  comment.  The  interaction  must  satisfy 
conservation  of  angular  momentum  and  energy 

(2 tiE, )xnbi  =  (2 tiEf  )',2b,  +  Jh /2ir,  ( 17) 

E,  +  EE=Ef  +  E,  +  BJ2,  (18) 

where  n  is  the  reduced  mass,  E,  and  Ef  are  the  initial  and 
final  center-of-mass  collision  energies,  b,  and  bf  are  the  ini¬ 
tial  and  final  impact  parameters,  EE  is  the  exoergicity  of  the 
process,  and  BJ 2  approximates  the  product  rotational  energy 
after  the  collision.  From  Eq.  (17),  large  values  of  J  can  be 
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achieved  through  highly  attractive  collisions  with  large  ini¬ 
tial  impact  parameters  and  small  final  impact  parameters 
with  the  concurrent  requirement  that  the  final  center-of- 
maas  kinetic  energy  is  small.4’  This  then  suggests  a  long- 
range  attraction  between  N2(a')  and  CO,  similar  to  that  in¬ 
ferred  for  N2(a)  +  CO.15  The  implication  from  surprisal 
analysis,  that  there  are  at  least  two  dynamic  mechanisms  for 
reaction  ( 14b),  further  suggests  a  strong  anisotropy  in  the 
potential  energy  surface,  by  analogy  to  observations  of 
CO(vJ)  formed  in  the  quenching  of  Na(3  2/>)  by  CO." 
Clearly,  detailed  potential  energy  surface  calculations  are 
required  to  test  this  hypothesis. 

CONCLUSIONS 

We  have  observed  extensive  rovibrational  excitation  in 
carbon  monoxide  from  quenching  of  discharged  nitrogen  in 
a  low  pressure  cryogenic  apparatus.  Analysis  of  the  data  has 
identified  a  contribution  to  CO(  t>  =  1 ,2 )  of  the  80  K.  compo¬ 
nent  from 

Nj(i>)  +  CO— N2(u  -  1 )  +  CO(t>  =  1 ) 

N2(i>)  +  CO-N2(u  —  2)  +  CO(u  =  2). 

Based  on  kinetic  and  energetic  arguments,  we  have  deter¬ 
mined  the  energy  transfer  process 

N2(a'  '2~  )  +  CO-N2(*,o)  +  COM  'll) 

-N2(*,t0  +  CCW) 

to  be  responsible  for  the  observed  rovibrational  excitation. 

The  energy  transfer  collisions  producing  the  rotational- 
ly  excited  component  exhibit  complex  dynamic  behavior, 
two  dynamic  mechanisms  have  been  identified,  one  respon¬ 
sible  for  low  vibrational  levels  and  another  that  contributes 
to  all  accessible  vibrational  levels.  Rigorous  potential  energy 
surface  calculations  will  be  required,  however,  to  determine 
the  true  dynamics  of  the  N2(fl‘)  +  CO  interaction. 
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The  ratios  of  overtone  and  fundamental  vibrational  Einstein  coefficients  for  NO(2T  'll )  have 
been  measured  by  spectrally  resolved  infrared  chemiluminescence  near  2.7-3.3  ftm  and  S.2-6.8 
ft m.  The  reactions  of  N  ( 2D,*P)  with  02,  in  the  presence  of  a  small  background  of  He  in  a 
cryogenic  low-pressure  reactor,  generated  vibrationally  excited,  rotationally  cold  (60  K) 

NO(d),  whose  emission  spectra  were  recorded  with  high  spectral  resolution.  Least-squares 
spectral  fitting  analysts  of  the  observed  overtone  and  fundamental  spectra  gave  vibrational 
bimd  intensities,  whose  ratios  at  each  emitting  vibrational  level  if  yielded  the 
(Au  rm  2)/(  Av  »  1)  Einstein  coefficient  ratios  for  i/  —  2-13.  The  results  provide  comparisons 
to  previous  theoretical  and  experimental  data,  and  reflect  the  behavior  of  the  dipole  moment 
function  for  NO(2f  3ri).  The  measured  ratios  indicate  an  overtone  Einstein  coefficient 
*  0-94  ±0.11  s  ” 1  for  an  assumed  fundamental  value  A 10  =  13.4  s” '. 


L  INTRODUCTION 

The  nitric  oxide  molecule  playa  key  roles  in  the  chemis¬ 
try  and  radiative  behavior  of  air-breathing  combustion  sys¬ 
tems,  chemical  lasers,  discharge  plasmas,  and  the  earth’s  at- 
In  all  of  these  applications,  the  infrared 
vibrational  transition  probabilities  for  NO  are  important  for 
diagnostic  applications  and  investigations  of  energy/radia¬ 
tion  transfer.  Owing  to  the  high  temperatures  and  nonequi- 
tibrium  conditions  in  which  NO  radiation  is  often  encoun¬ 
tered,  the  emitting  species  can  be  highly  internally  excited, 
requiring  knowledge  of  vibration-rotation  transition  mo¬ 
ments  at  large  internuclear  separations.  This  information  for 
NO  is  available  only  from  theory,  with  little  or  no  experi¬ 
mental  validation.  In  this  paper,  we  present  experimental 
data  for  the  branching  ratios  of  fundamental  and  first  over¬ 
tone  emission  from  highly  vibrationally  excited  NO,  as  ob¬ 
served  in  a  cryogenic  chemiluminescence  reactor.  These 
data  provide  insight  into  the  shape  of  the  electric  dipole  mo¬ 
ment  function,  as  well  as  an  indirect  determination  of  the 
absolute  (o’  =  2)  —  (v*  —  0)  overtone  transition  probabili¬ 
ty. 

In  the  specification  of  absolute  emission  intensities  and 
spectral  distributions,  one  requires  the  vibrational  emission 
rates  (Einstein  coefficients,  A^ ),  averaged  over  all  ther¬ 
mally  populated  rotational  levels,  for  all  the  vibrational 
states  of  interest.  While  many  ground  state  absorption  mea¬ 
surements  have  determined  values  for  the  (v',vm)  =  (1,0) 
and  (2,0)  transition  strengths,  the  scaling  of  Ap.t.  with  v'  is 
only  partially  understood.  Previous  empirical  determina¬ 
tions’-3  of  the  dipole  moment  near  equilibrium  internuclear 
separation  provided  limited  information  for  low  v'.  In  a  de¬ 
finitive  theoretical  treatment,  Billingsley3”6  performed  an  ab 
initio,  multiconfiguration  self-consistent-field  (MCSCF) 
calculation  of  the  dipole  moment  function  of  NO(>f  3II), 
and  rigorously  evaluated  the  rotationally  dependent  and 


thermally  averaged  Einstein  coefficients  for  the  fundamen¬ 
tal  (Auk  1)  and  first  overtone  (Av  =*  2)  transitions  over  a 
large  range  of  v' and./'  (if  =  1-20^'  -  0.5-33.5).  The  over¬ 
lap  of  the  ah  initio  dipole  moment  function  and  the 
NO(*2Il)  potential  is  illustrated  in  Fig.  1.  The  absolute 
values  and  relative  scalings  given  by  this  work  have  been 
widely  utilized;  however  the  absolute  values  of  are  now 

known  to  be  erroneously  small.7 

In  the  only  previous  experimental  investigation  of  vibra¬ 
tional  dependence.  Green  era/.*  determined  Einstein  coeffi¬ 
cient  branching  ratios  for  NO(v'  =  2-9)  from  Av  =  1  and 
Au  =  2  fluorescence  spectra  in  electron-bombarded  N2/02 
mixtures.  These  measurements  confirmed  Billingsley’s6  pre¬ 
dicted  vibrational  scaling  of  the  branching  ratios  for  if  =*  3- 
7,  and  established  that  the  dipole  moment  function  of  Mi¬ 
chels1  was  incorrect.  However,  the  experimental  data  leave 
considerable  uncertainty  for  the  (2,0)/(2,l )  ratio,  and  are 
not  definitive  for  if>8.  Nevertheless,  the  data  of  Green  et 
al.*  verify  the  shape  (second  derivative)  of  the  ah  initio  di¬ 
pole  moment  function  at  intermediate  nuclear  separation. 
Thus,  a  reasonable  approach  for  selecting  values  is  to 

apply  the  Billingsley6  vibrational  scalings  to  absolute  values 
ofAuo  and/or  A M  determined  from  careful  absorption  mea¬ 
surements. 2’9,10  This  approach  was  adopted  by  Rothman  et 
al.1  to  evaluate  absorption  strengths  for  both  the  fundamen¬ 
tal  and  the  first  overtone;  the  resulting  absolute  values  were 
slightly  modified  by  Rawlins  et  al."  for  the  analysis  of 
NO(  Av  =  1 )  chemiluminescence  spectra. 

We  have  previously  reported"  detailed  laboratory  in¬ 
vestigations  of  the  vibrational  chemiluminescence  from  nas¬ 
cent  NO(ivf)  formed  by  reactions  of  metastable  atomic  ni¬ 
trogen  with  O, ,  as  originally  observed  by  Kennealy  et  al. 12 
In  the  recent  experiments,1 1  extensive  rotational  and  vibra¬ 
tional  excitation  was  observed  in  nitric  oxide  formed  under 
nearly  collisionless  conditions  by  the  reactions 
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N(*D)  +  O,  -NO(»)  +  O  0D*P),  (RI) 

N(*l*)  +  O,  -NO(o^)  +  0(  'S^D’P).  (R2) 

The  NO(v)  distribution*  from  (Rl)  peak  near  u  =  7,  and 
are  detectable  to  o  «  14.  Reaction  (Rl)  ia  responsible  for 
auroral  cbctniexchation  of  NO^),'*14  and  for  most  of  the 
production  of  NO  in  the  earth’s  thermosphere. 13 

Reaction  (R2)  produces  NO(tv/)  with  a  high  degree  of 
rotational  excitation  in  each  of  several  (at  least  eight)  vibra¬ 
tional  levels.11  The  spectral  data  exhibit  distinctive  X- 
branch  band  heads,  which  signify  substantial  populations  of 
rotational  states  in  the  range/ ~  60.5-120.5  (i.e.,  >leVof 
rotational  energy).  Similar  band  heads  have  been  observed 
in  high-altitude  auroral  spectra  ofNO  emission, 14  indicating 
a  significant  role  for  (R2)  and/or  similar  chemiexcitation 
processes  in  the  particle-bombarded  upper  atmosphere. 

In  the  analysis  of  spectral  distributions  for  the  high  vi¬ 
brational  and  rotational  states  populated  by  reactions  (Rl) 
and  (R2),  it  is  necessary  to  probe  the  dipole  moment  func¬ 
tion  over  a  much  larger  range  of  the  NO  potential  surface 
than  has  previously  been  examined  experimentally.  In  addi¬ 
tion,  it  is  often  important  to  specify  the  overtone  spectral 
distribution  corresponding  to  that  of  the  fundamental,  or 
vice  versa.  To  these  ends,  we  have  used  the  infrared  chemilu¬ 
minescence  approach,  previously  applied  to  determine  tbe 
nascent  product  distribution  for  reactions  (Rl)  and  (R2), 
to  simultaneously  observe  At;  —  land  Acs  2  spectral  distri¬ 
butions  from  those  reactions.  To  eliminate  complications 
due  to  multimodal  rotational  distributions  and  spectral  band 
overlap,  we  introduced  low  levels  of  helium  bath  gas  which 
caused  extensive  rotational  cooling  and  slight  vibrational  de¬ 
activation  of  the  initial  state  distributions.  This  effect,  com¬ 
bined  with  the  cryogenic  temperature  of  the  reaction 
chamber,  resulted  in  extremely  low-temperature  Boltzmann 
rotational  distributions,  which  permitted  accurate  detenni- 
nstionsofthe  (Ac  ~  2)/(Ac  »  1)  branching  ratios  over  the 


range  c'  —  2-12.  The  results  clarify  previous  determinations 
at  low  v\  and  indicate  significant  departure  from  Billings¬ 
ley’s*  scaling  at  high  o’. 

A  subsequent  paper17  will  address  the  analysis  of  the 
branching  ratios  to  determine  an  empirical  dipole  moment 
function  for  NO  and  the  corresponding  Einstein  coefficients. 
We  describe  here  the  experimental  methods,  spectral  analy¬ 
sis,  and  branching  ratio  results  as  compared  with  previous 
work. 

IL  EXPERIMENTAL  MEASUREMENTS 

The  experiments  were  performed  in  the  cryogenic  CO¬ 
CHISE  (cold  chemiluminescent  infrared  stimulation  exper¬ 
iment)  facility  at  the  Geophysics  Directorate.  The  design 
and  operation  of  this  facility  are  described  in  detail  else¬ 
where;1*  the  measurement  conditions  were  essentially  the 
same  as  those  for  the  previous  NO  chemiluminescence  inves¬ 
tigation.1’  In  brief,  the  reaction  cell  and  surrounding  radia¬ 
tive  environment  are  cooled  to  approximately  15-20  K, 
which  eliminates  background  radiation  and  provides  rapid 
cryopumpiiig  of  the  reagent  gases.  The  infrared  detection 
system  consists  of  a  cryogenic,  scanning  grating  monochro¬ 
mator  and  a  hquid-behum-cooled  Si:As  detector,  focused  to 
infinity  along  the  axis  of  the  reaction  cell.  Reagent  gases,  at 
80  K,  enter  the  cell  through  a  aeries  of  four  opposing  jets,  and 
mix  along  the  cell  axis  as  iDustrated  in  Fig.  2.  In  the  present 
case,  nitrogen  metastabies,  generated  by  microwave  dis¬ 
charge  of  flowing  N,  /At  (12%  N,)  mixtures  near  1  Tour, 
enter  from  one  side  of  the  cell  and  mix  along  the  centerline 
with  a  mass-balanced  opposing  flow  of  Of  This  results  in  a 
rapid  chemical  reaction  producing  chemiexcited  NO (vj), 
which  radiates  in  infrared  bands  near  5.4  (Av  *  1)  and  2.7 
(turns!)  fim.  Tbe  discharges  are  modulated  with  a  23  Hz 
square  wave,  and  the  emission  from  the  reaction  zone  is  ob¬ 
served  via  phase-sensitive  detection. 

In  other  applications  of  this  technique,1*'21  the  pressure 
in  the  reaction  zone  is  typically  3-5  mTorr,  approaching 
single  collision  conditions.  For  the  present  experiments,  be- 


o* 


HO.  2.  Diagram  of  COCHISE  reaction  chamber. 
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Bum  wm  admitted  to  achieve  a  steady  state  (static)  pressure 
of  10  mTorr,  quenching  out  the  high-/  Mates.  Preliminary 
measurements  in  the  absence  of  He  gave  overtone  spectra 
which  were  difficult  to  analyze  unambiguously,  owing  to  the  | 
extensive  spectral  overlap  between  poorly  resolved  high-/  ^ 
and  thennalized  rotational  components  of  adjacent  vibra-  T 
tional  bands.  The  addition  of  small  amounts  of  He  resulted  | 

in  significant  rotational  cooling,  and  a  pressure  of  10  mTorr 
was  sufficient  to  eliminate  the  high-/  components.  This  is 
demonstrated  by  the  absence  of  /i -branch  band  heads  in  the 
Av  =  1  spectra,  and  by  the  conformity  of  the  rotational  dis¬ 
tributions  to  a  single-temperature  Boltzmann  form  for  all 
vibrational  levels. 

Spectra  of  the  Au  -  1  and  Av  «  2  emission  were  record¬ 
ed  sequentially,  as  matched  pain,  using  order-sorting  inter¬ 
ference  filters  to  isolate  each  band.  A  total  of  four  matched 
pairs  were  recorded,  all  for  essentially  the  same  conditions  of 
pressure,  temperature,  flow  rates,  and  discharge  operation. 

All  spectra  were  obtained  at  a  resolution  of 0.0067  fim  (foil 
width  at  half-maximum),  corresponding  to  approximately 
2.3  cm" 1  for  the  fundamental  band  and  9.2  cm  * 1  for  the 
overtone  band.  The  uncertainty  in  the  observed  wavelengths 
is  ±  0.003  fim  due  to  a  periodic  fluctuation  in  the  mono¬ 
chromator  scan  rate  at  cryogenic  temperature. 

Before  and  after  the  spectral  measurements,  the  optical 
system  was  calibrated  for  absolute  and  relative  spectral  re¬ 
sponse  using  a  black  body  radiation  source  imbedded  in  the 
mirror  at  the  end  of  the  reaction  cell.  The  calibrations  em¬ 
ployed  btackbody  temperatures  between  350  and  400  K, 
with  a  temperature  measurement  uncertainty  of  ±  3  1C 
This  results  in  a  systematic  uncertainty  of  ±  12%  in  the 
measured  (Av>2)/(Av>e  1)  intensity  ratios. 

A  representative  spectral  pair  is  shown  in  Fig.  3,  with 
the  NO  transition  band  centers  labeled.  Note  the  logarith¬ 
mic  intensity  scale  illustrating  a  dynamic  range  of  approxi¬ 
mately  two  orders  of  magnitude.  The  NO(At; »  1)  system  is 
clearly  resolved  up  to  rf  «  13,  with  no  apparent  interfering 
radiators.  However,  several  additional  features  appear  with 
the  NO(Ac  —  2)  system.  The  (5,2)  and  (2,0)  bands  of  the 
N*  (  W*tka  —  B  JIlg )  electronic  transition22  are  evident  near 
2.6  and  3.3  ftm,  respectively,  the  3.3  pn  band  ultimately 
limits  the  range  of  for  which  the  analysis  can  be  carried 
out  In  addition,  several  Rydberg  transitions  of  Ar  are  ob¬ 
servable  as  scattered  light  from  the  discharges. u  In  particu¬ 
lar,  a  weak  Ar  fine  near  2.69  fim  partially  obscures  the 
NO(2—0)  emission,  necessitating  a  factor  of  s:  2  downward 
correction  in  the  observed  intensity  of  this  band. 


111.  ANALY8I8  ANO  RESULTS 

The  spectra  were  analyzed  using  a  linear  least-squares 
spectral  fitting  method  which  we  have  employed  extensively 
in  the  past. 1,1  The  NO  line  positions  and  strengths 
were  calculated  with  the  spectroscopic  constants  given  in 
Rawlins  et  of."  (as  modified  slightly  from  those  of  Gold¬ 
man  and  Schmidt24 )  using  the  assumption  of  Honl-London 
scaling  of  the  rotational  line  strengths.23  From  the  predic¬ 
tions  of  Billingsley,3’*  the  Hdnl-London  assumption  should 
be  accurate  for  the  low-/  levels  sampled  in  these  measure- 
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ments.  In  particular,  spin-uncoupling  and  vibration-rota¬ 
tion  interactions  will  not  impact  the  analysis  of  the  data. 
Matching  of  the  observed  rotational  envelopes  to  the  com¬ 
puted  spectra  gives  a  uniform  rotational  temperature  of  60 
K,  signifying  that  the  presence  of  He  causes  significant  ther¬ 
mal  coupling  between  the  80  K  reagent  gases  and  the  20  K 
reactor  wall.  At  60  K,  the  third  rotational  level  carries  the 
maximum  population,  and  99%  of  the  population  lies  in  the 
lowest  ten  rotational  levels.  Examples  of  the  computed  spec¬ 
tral  fits  are  shown  in  Figs.  4  and  5. 

The  least-squares  solutions  to  the  spectral  fitting  proce¬ 
dure  are  the  integrated  intensities  for  each  vibrational  transi¬ 
tion,  [NO(u')  as  (dotted  in  Fig.  6  for  a  matched  pair 
of  overtone  and  fundamental  spectra.  By  taking  ratios  of 
these  values  at  each  &  for  a  given  spectral  pair,  we  determine 
the  Einstein  coefficient  ratios  Ayy  _  t/A,.y  _ , ,  which  rep¬ 
resent  the  branching  ratiosbetween  Ac  =  2and  Au  =  1  tran¬ 
sitions.  The  observed  Einstein  coefficients  are  thermally 
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FIG  4  Data  ( light  line,  shown  with  best  It 
and  separately  above)  aad  least -squares  St 
(dark  taw)  to  NO(Ag»  PemiaooB  spec¬ 
trum.  Rotatsooai  tcaipcnture  aad  spectral 
mol  ik  ion  art  60  K  aad  0.0067  ftm.  The 
standard  deviation  of  the  spectral  It  is 
0041. 


avenged  at  60  K.  However,  since  Hdnl-London  scaling  ap¬ 
plies  (no  spin  uncoupling),  these  Einstein  coefficients  an 
equivalent  to  the  rotationless  Einstein  coefficients  of  the 
molecule,  and  an  independent  u  temperature  to  at  least  <00 
K  as  shown  by  Billingsley.*  The  observed  ratios  for  ail  the 


FKI.  S.  Data  (Isht  hat  shows  with  the  beat  St  aad  separately  above)  aad 
least -aqaaraett  (dart  has)  to  NO<  As-  2)  eansaion  ipactrvai  Rotational 
fnpantnrv  aad  spactrsl  rseotathm  are  <0  K  aad  0.0067  jam.  The  slight 
nwnMchet  aaar  2.9  aad  3. 1  /am  malt  from  aontawarify  fet  the  rrynptlc 
(raaias  drive.  The  staadard  deviation  at  the  spectral  St  is  0.1 5. 


data,  and  avenge  values  for  each  o',  an  tabulated  in  Table  I. 
The  uncertainties  given  in  Table  I  an  statistical  errors  deter- 
mined  for  the  fitting  and  avenging  procedures,  and  do  not 
incorporate  the  systematic  uncertainty  of  the  (dative  re¬ 
sponse  calibration. 


Upper  VSaaSoiul  Level  (V) 


FIQ.  6.  Vibratioaal  band  inleasitics  determined  front  a  least-squares  6t  of 
thadaai ratal  anil  overtone  spectra.  Error  bars  denote  ±  la  standard  devi¬ 
ation. 
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TABLE  L  ( At  -  2)/<  A»  «  1 )  bwd»q  nbo*. 


Weighted 

4  A,y  .  a/^v  -  I  <  ±  •»)  avenge  (  ±  l<r) 


2 

oo4i±ooi7 
0.031 ±0026 
0.034  ±0.011 
0.036  ±  0.016 

0037  ±0001 

3 

aoto±aou 

0.014  ±  0016 
0.096  ±0.009 
0.096  ±  0.020 

0079  ±0.007 

4 

0.096  ±  0.007 
ai07±aaos 
0.107  ±0.003 

ao96±aoio 

0.103  ±0003 

5 

a  131  ±0.006 
0120  ±0006 
0.152  ±0.003 
at3i±aoo9 

0.136  ±0.009 

6 

0.152  ±0.006 
ai46±aoo6 
0.172  ±0.004 
ai47±aooi 

0159  ±0.007 

7 

0.193  ±0.003 
0.215  ±0.010 
0.221  ±0.007 
0.163  ±0.010 

0.202  ±0.012 

S 

0.202  ±0lOU 
0215  ±0.013 
0247  ±0.012 
0190  ±0016 

0.214  ±0.014 

9 

0293  ±0034 
0331  ±0042 
0339  ±0.031 
0263  ±0036 

0313  ±0.022 

to 

0.2*9  ±0.066 
0239  ±0063 
0.379  ±0.074 
0.263  ±0073 

0293  ±0027 

11 

0390  ±0234 
03!9±0197 
0439  ±0.222 
0423  ±0.263 

0394  ±0.029 

12 

0397  ±0.357 
0313  ±0333 
0279  ±0204 
0223  ±0.267 

0.296  ±0.030 

13 

0572  ±0.732 
0299  ±  1.001 
0319  ±0.343 
0092  ±0.406 

0262  ±0.099 

IV.  DISCUSSION 

A.  Comparison*  to  previous  dita 

The  equations  relating  Einstein  emission  coefficients  to 
the  dipole  moment  function  are  discussed  in  detail  by  Bil¬ 
lingsley*  The  band  Einstein  coefficients  A  i,.r,  represent 
thermal  averages,  for  a  given  rotational  temperature,  of  the 
individual  state-specific  Einstein  coefficients  for  a  given  vi¬ 
brational  transition.  These  state-specific  Einstein  coeffi¬ 
cients  are  proportional  to  the  square  of  the  dipole  moment 
integral  connecting  the  two  states, 


|  j£.  dl'(rV<(r)p*(r)dr|  , 

where  p'(r)  and  p*  (r)  are  the  wave  functions  of  the  upper 
and  lower  states  andp(r)  is  the  dipole  moment  as  a  function 
of  intern uclear  distance  r.  For  an  electrically  harmonic  sys¬ 
tem,  wherep(r)  is  linear,  the  magnitudes  of  the  At;  =  1  tran¬ 
sition  probabilities  are  determined  from  the  first  derivative 
of  the  dipole  moment  function,  while  those  for  the  overtone 
transitions  are  determined  solely  by  the  mechanical  anhar- 
monicity  of  the  system,  i.e.,  the  nonorthogonality  of  the 
wave  functions  p*  and  p*.  As  the  dipole  moment  function 
acquires  curvature,  increases  in  the  second  and  third  deriva¬ 
tives  strongly  increase  the  overtone  transition  probabilities. 
In  the  case  of  OHCJf2Il),  the  dipole  moment  function  is 
strongly  curved  near  equilibrium  in  tern  uclear  separation, 
and  as  a  result  the  Av  =  2  Einstein  coefficients  exceed  those 
for  Av  *  1  at  unusually  low  i/.2*00 

For  NO(jrJI1),  the  dipole  moment  function  is  small 
and  only  mildly  curved  near  equilibrium  separation.2'4  Thus 
we  expect  the  Ao  =  1  Einstein  coefficients  to  be  primarily 
sensitive  to  the  slope  of  fifr),  and  fairly  insensitive  to  its 
shape.  The  scaling  of  ^  v  _ ,  with  o'  should  be  close  to  har¬ 
monic  (i.e.,  a.  o')  near  the  bottom  of  the  potential  well,  and 
will  tail  off  at  higher  o'  due  to  the  increasing  effects  of  me¬ 
chanical  anharmonicity  in  the  wave  functions.  In  addition, 
increasing  curvature  in  ft(r)  at  large  r  will  affect  the  Av  —  1 
values  at  some  high  (but  unknown)  i/.  In  contrast,  even 
modest  curvature  in  /*  ( r )  will  profoundly  affect  the  absolute 
values  and  o‘ '  scaling  of  the  Av  =  2  Einstein  coefficients. 
Thus  the  measurement  of  (Av  m  2)/(Av  =  1)  branching  ra¬ 
tios  over  a  wide  range  of  o'  provides  information  on  the  cur¬ 
vature  of  the  dipole  movement  function  over  a  large  range  of 
r. 

The  average  ( Av  =  2)/(  Av  —  1 )  ratios  are  {dotted  in 
Fig.  7,  together  with  previously  measured  and  predicted  val¬ 
ues.  The  agreement  with  Green  et  al.*  is  excellent  over  the 
range  4<v'<7,  where  their  data  have  reasonable  precision. 
The  COCHISE  measured  values  deviate  significantly  from 
the  Billingsley  scaling  at  both  low  and  high  o',  and  indicate 
the  presence  of  a  broad  maximum  near  o'  =  11.  The  exis¬ 
tence  of  such  a  maximum  could  have  significant  implica¬ 
tions  for  the  shape  of  the  dipole  moment  function  at  large 
intern  uclear  separation  and  for  the  behavior  of  at  high 
o'  and/or  Unfortunately,  due  to  the  limited  precision 
(low  signal/noise)  of  our  spectra  in  this  region,  we  cannot 
obtain  sufficiently  accurate  o'  =  13  data  to  define  the  appar¬ 
ent  downturn  in  the  ratios  for  v>  12.  Nevertheless,  the  data 
clearly  show  significant  deviations  from  the  shape  of  the  ab 
initio  dipole  moment  function  for  both  equilibrium  and  large 
intemuclear  separations. 

The  ratio  values  in  Table  I  are  appropriate  for  rotational 
temperatures  below  =;  100  K,  and  are  probably  reliable  for 
temperatures  up  to  at  least  300  K.  Billingsley’  *  found  that 
the  individual  line  strengths  departed  only  slightly  from 
Hbnl-London  scaling,  and  as  a  result  the  thermally  aver¬ 
aged  Einstein  coefficients  were  virtually  independent  of  tem¬ 
perature  up  to  600  K.  The  extension  of  these  values  to  higher 
J‘  and  temperature  requires  further  investigation.  At  some 
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FIO.  7.  Cowpirttaa  of  cxperimaHal  branching  ratio*  to  (»)  prtviou*  me»- 
Mmnwwund  (b)  prediction*  from  dipole  moment  fiwctioM.  Solid  circk*: 
we^Med  averaga  froai  Table  I,  with  error  ban  tijnifying  ±  ltriundard 
deviation  from  the  weighted  mean.  Open  circle*:  data  from  Green  «  a/.. 
Ref.  t.  Solid  carve:  BiKagriey,  Ref.  6.  Deshed  curve:  Chandraiah  and  Cho, 
Ref.  2  (calculated  in  Ref.  6).  Dotted  curve:  Michels,  Ref  I  (calculated  in 
Ref.  6). 


temperature,  will  become  sensitive  to  rotational  distri¬ 

bution,  since  the  scaling  of  A  (/')  diverges  from  Honl-Lon- 
doo  as  /’  increases.3  It  is  possible  that  the  high-/  distribu¬ 
tions  resulting  in  R -branch  band  heads"  are  significantly 
affected,  since  the  effects  of  spin-uncoupling  and  vibration- 
rotation  interactions  may  be  substantial,  and  further  since 


these  states  sample  the  extremes  of  internudear  separation 
where  curvature  in  the  dipole  moment  function  is  more  pro¬ 
nounced.  Further  analysis  of  the  COCHISE  data  to  infer  an 
empirical  dipole  moment  function,  examine  high-/  effects, 
and  compute  Au,c.  ( T)  will  be  reported  elsewhere.17 

8.  Abuohrtn  Einstein  cosflldsnte:*,*  and4M 

Most  of  the  experimental  data  base  on  NO(  v)  transition 
probabilities  comes  from  absorption  measurements  on  the 
(1,0)  and  (2,0)  transitions.  A  list  of  previously  published 
experimental  and  theoretical  values  for  the  absorption 
strengths  and  Einstein  coefficients  for  these  transitions  is 
given  in  Table  II.  The  measurements  for  the  ( 1,0)  transition 
fall  into  three  groups:  those  with  Aija  =  13-14  s~ ",  those 
with  12-13,  and  a  few  early  measurements  of  7-8.  The  work 
of  Holland  et  a/.,10  together  with  the  extrapolated  line-by¬ 
line  data  of  Mandin  et  al?  as  reported  by  Rothman  et  al.,1 
give  consistent  values  of4l>0  of  13.6  and  13.3  s'1,  respective¬ 
ly.  A  discussion  comparing  these  values  to  previous  mea¬ 
surements  can  be  found  in  Holland  etal. 10  In  brief,  the  high 
spectral  resolution  and  higher  measurement  precision  of 
Reft.  9  and  10  indicate  a  significant  discrepancy  with  the 
previous  lower  values,  and  strongly  support  a  value  for  Aixt 
between  13  and  14s~'.  Indeed,  Rothman  etal.1  invoked  the 
extrapolated  line-by-line  data  of  Mandin  el  al.9  in  their  com¬ 
pilation  of  atmospheric  line  absorption  parameters.  We  have 
adapted  the  value  AlJ0  =  13.4  s~ 1  at  300  K.  This  value  is 
further  supported  by  recent  shock  tube  measurements  of 
spectrally  resolved  NO(Av  ™  1)  emission  between  800  and 
2500  K,  which  give .4I  0  *  13.2 s~'  (  ±  10%).2' 

The  (2,0)  transition  probability  is  less  well  determined. 
Experimental  values  range  from  0.8  to  1.1  s~ '.  The  most 
careful  measurements  are  those  of  Chandraiah  and  Cho,2 
whose  empirical  dipole  moment  function  provides  good 
agreement  with  the  COCHISE  branching  ratios  at  low  if 
(see  Fig.  7).  We  note,  however,  that  their  result  for  the  ( 1,0) 
band  appears  to  be  about  10-12%  low. 

The  COCHISE  data  give  an  accurate  value  of  A2J0/A^ 

(  ±  2<r  =  7.9%, cf.TabieI),whichcanbeusedtodetermine 
A2M  with  good  precision.  Billingsley6  gives  the  ratio 
—  190, 1.89,  and  1.92  for  the  dipole  moment  func¬ 
tions  of  Reft.  4, 2,  and  1,  respectively.  Since  the  dipole  mo¬ 
ment  function  is  nearly  linear  near  equilibrium  internudear 
separation,  this  particular  ratio  is  dominated  by  the  mechan¬ 
ical  anharmonidty  of  NO  (i.e.,  the  radial  wave  functions  for 
o  =  0,  1,  2),  and  is  relatively  independent  of  the  details  of 
slope  and  shape  of  the  dipole  moment  function.  Thus  the 
,4 2,1  /4  i,o  ratio  is  s  1.9  for  three  very  different  dipole  mo¬ 
ment  functions,  and  is  probably  quite  accurate.  We  can  then 
determine  from  the  experimental  data  the  ratio 
AM/AlJ0  =  0.070,  which  is  identical  to  that  determined  by 
Chandraiah  and  Cho2  and  about  1.6  times  that  computed  by 
Billingsley.6  {In  comparison,  the  experimental  value  for  this 
ratio  is  2.4  times  larger  than  for  the  electrically  harmonic 
(linear  pir)  ]  case.}  Application  of  our  recommended  value 
for  Alj0  gives  the  determination 

A2j0  =0.94  ±0.11  s', 

where  the  indicated  error  is  the  systematic  uncertainty  in  the 
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TABLE  11. 


ikti —  1  and  -  2  transitions  of  NO. 


5  6  7  8  9  10  11  12  13 

UpparVVxabonai  Laval  (v") 


FIG.  8.  Relative  scaling  of  Au  —  I  Einstein  coefficients  with  p'.  Solid  curve: 
Billingsley,  Ref.  6.  Dashed  curve;  Chandraiah  and  Cho.  Ref.  2  (calculated 
in  Ref.  6).  Dotted  curve:  Michels,  Ref.  1  (calculated  Hi  Ref.  6).  Dash/dot 
curve;  computed  for  a  linear  dipole  moment  function. 
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relative  response  calibration.  This  value  is  12%  higher  than 
(but  within  experimental  uncertainty  of)  that  measured  by 
Chandraiah  and  Cho,2  and  is  a  factor  of  2  larger  than  pre¬ 
dicted  from  the  dipole  moment  function  of  Billingsley.*  The 
relatively  large  value  for  ^2.0  indicates  that  the  dipole  mo¬ 
ment  function  has  mote  curvature  near  equilibrium  intemu- 
clear  separation  than  was  predicted  from  the  oh  initio  theory. 

Potentially,  the  assumption  of  Billingsley  scaling  of 
Ae.y  _ ,  can  be  used  with  the  COCHISE  data  to  infer 
Ajy  _  2  values  to  higher  t»';  however,  this  procedure  will 
have  increasing  uncertainty  at  higher  o'  as  curvature  in  the 
dipole  moment  function  plays  more  of  a  role  in  the  o'  depend¬ 
ence  of  A'.y  _ , .  To  evaluate  the  validity  of  this  procedure, 
we  examine  the  predicted  Aa.  u.  _ ,  scalings  as  shown  in  Fig. 
8.  Predictions  for  the  dipole  moment  functions  of  Refs.  1, 2, 
and  4  were  taken  from  Ref.  6.  In  addition,  we  have  plotted 
the  predicted  scaling  for  a  linear  dipole  moment  function  as 
computed  using  Rice-Ramsperger-Kassel  (RKR)  wave 
functions  for  NO(T2n).  (While  the  absolute  values  of 
A".  *.  _ ,  depend  on  the  assumed  slope  of  ft(r),  the  relative 
values  Av.  v.  _l/At  0  do  not.)  It  is  clear  from  Fig.  7  that  the 
relative  A  factors  for  the  fundamental  band  are  fairly  insensi¬ 
tive  to  modest  variations  in  ^(r)  curvature  for  u'<4.  Fur¬ 
thermore,  the  close  agreement  between  our  A20/Alja  ratio 
and  that  of  Chandraiah  and  Cho,2  together  with  the  similar¬ 
ity  in  o'  scaling  of  the  A  factors  for  to  —  1  computed  by 
Billingsley6  from  the  ab  initio*  and  empirical2  dipole  mo¬ 
ment  functions,  suggest  that  the  Billingsley  scaling  is  indeed 
reliable,  at  least  up  to  o'  —  6. 

With  the  caveat  that  curvature  in  jj(r)  at  large  r  can 
affect  the  >4^  _,/d,  0  scaling  at  high  0,  we  tentatively  con¬ 
clude  that  Billingsley  scaling  is  a  good  approximation  for  the 
fundamental  band,  with  an  uncertainty  increasing  with  o’  to 
perhaps  ±  20%  near  v'  —  12.  (However,  we  note  that  there 
is  no  firm  basis  for  extending  this  scaling  to  higher  o',  say 
o'  ~  20. )  Accordingly,  we  have  applied  this  scaling  of 
A".y  _ ,  to  estimate  values  forA^.y  _  2  using  a  quartic  least- 
squares  fit  to  the  branching  ratio  data  of  Table  I.  The  results 
are  luted  in  Table  III.  Given  the  combined  uncertainties  in 


TABLE  III.  Einstein  confidents  for  NO(Ap»  2)  transitions. 


d 

Estimated  from 
COCHISE 

Ab  initio 
(Ref.  6) 

2 

0.94  (  ±  12%) 

0.46 

3 

2.74 

1.51 

4 

4.8 

3.10 

5 

7.1 

4.90 

6 

9.8 

7.29 

7 

13.0 

9.63 

8 

16.6 

12.50 

9 

20.2 

15.68 

10 

23.9 

19.14 

II 

26.5 

22.94 

12 

27.9 

27.15 

13 

27.2  (±40%) 

31.88 

the  branching  ratio  measurements  and  the  dt,  (  scaling, 
we  conservatively  estimate  uncertainties  in  Av,  a.  _  2  ranging 
from  ±  12%  for  o'  =  2  to  ±  40%  for  o'  =  13. 

V.  SUMMARY  AND  CONCLUSIONS 

We  have  used  the  infrared  chemiluminescence  tech¬ 
nique,  in  a  cryogenic  reaction  vessel,  to  measure  the  over¬ 
tone/fundamental  branching  ratios  for  vibrational  emission 
from  NO(2f2n,i>).  A  combination  of  low  rotational  tem¬ 
perature  and  high  spectral  resolution  eliminated  complica¬ 
tions  of  non-Boltzmann  rotational  distributions,  and  mini¬ 
mized  spectral  overlap  between  adjacent  bands.  The  spectral 
resolution  also  reduced  complications  in  the  analysis  due  to 
interference  from  other  radiating  species.  The  measure¬ 
ments  cover  ranges  of  2-13  in  o'  and  0-12  in  o',  and  extend 
the  previous  experimental  data  base  to  both  lower  and  high¬ 
er  vibrational  levels.  The  agreement  with  predictions  from  a 
previous  ab  initio  calculation  is  generally  favorable,  but  de¬ 
tailed  comparison  indicates  significant  differences  in  slope 
and  curvature  of  the  dipole  moment  function  at  equilibrium 
and  large  interaudear  separations.  These  differences  will  af¬ 
fect  the  magnitudes  and  o'  scaling  of  the  Au  =  1  and  Atr  =  2 
Einstein  coefficients,  the  more  significant  effect  being  on  the 
overtone  transitions.  An  experimental  value  derived  for  A2  0 
agrees  with  previous  determinations  and  shows  the  theoreti¬ 
cal  value  to  be  a  factor  of  2  too  low.  Further  details  of  these 
effects,  as  well  as  the  scaling  of  transition  strengths  for  high 
vibrational  and  rotational  excitation,  must  be  addressed 
through  derivation  of  an  empirical  dipole  moment  function. 
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A  definitive  measurement  of  the  rate  coefficient  for  the  quenching  of  N(l2))  by  Of1/*)  is 
reported.  The  0(  iP)  atoms  were  generated  by  titrating  NO  directly  into  the  active  nitrogen 
flow.  Analysis  of  the  results  required  that  the  rate  coefficient  for  the  reaction 
N(Ji))  +  NO— Nj  +  O  be  known  accurately,  and  this  was  also  determined.  A  finite  mixing 
time  correction  is  also  necessary.  The  best  estimate  of  the  rate  coefficient  from  this  work  is 
(6.9  +  0.7,  —  1. 1 )  x  10 '  *’  cm*  s  ~ 1  at  T  —  298  K,  considerably  smaller  than  a  previous 
measurement  [J.  Phys.  Chem.  92,  5977  ( 1988) ]  and  in  much  better  accord  with  values 
required  by  atmospheric  models. 


I.  INTRODUCTION 

N( ’/))  is  an  important  constituent  of  the  mesophere 
and  thermosphere.  Atmospheric  modeling  studies'  *  indi¬ 
cate  that  the  reaction 

a., 

N(2f»  +  Oj  — NO(i»)  +  O('P.'D)  (1) 

is  the  primary  source  of  the  NO  concentrations  first  ob¬ 
served  in  rockel  experiments.4  While  the  value  of  the  rate 
coefficient  for  reaction  ( 1 )  is  well  established,5',u  the  coeffi¬ 
cient  for  the  quenching  reaction 
A. 

N(-7»  +  0(7>)-N(4S)  +OCP)  (2) 

is  still  the  subject  of  some  dispute.  The  relative  magnitude  of 
these  two  rate  coefficients  is  important  in  determining  the 
partitioning  between  the  two  product  channels.  Modeling 
calculations  require  that  reaction  ( I )  be  the  dominant  loss 
channel  for  N (*/)).  A  recent  experimental  measurement" 
of  the  rate  coefficient  for  reaction  (2)  reported  a  value  of 
3.4x  10  "xexp(  —  145 /T)  cm’s" ',  which  would  corre¬ 
spond  to  2.1  X  10 " 11  cm’  s'  'at  T  —  300  K;  the  quoted  error 
limits  are  on  the  order  of  ±  50%.  This  value  is  considerably 
higher  than  previous  determinations  and  more  discordant 
with  preferred  modeling  values211'15  which  require 
ktt  <  10  cm’s"  If  this  rate  were  correct,  then  N(2£>) 
would  contribute  to  the  overall  depletion  of  NO  by  enhanc¬ 
ing  the  N(4S)  concentration  which  reacts  rapidly  with  NO 
through  the  reaction 

N(4S)  +NO-N2  +  0(’F).  (3) 

This  would  necessitate  a  re-evaluation  of  the  role  of  N(20) 
in  the  upper  atmosphere  and  the  identification  of  another 
source  term  for  NO. 

We  have  remeasured  the  rate  coefficients  for  both  reac¬ 
tions  ( 1 )  and  (2)  in  order  to  establish  the  true  partitioning 
between  these  two  N(*D)  loss  channels.  This  was  done  using 
the  FACEL1F  reactor"’  with  multiplioton  ionization 
(MPI)  detection  of  reactant  and  product  species.  The  ex¬ 


perimental  method  makes  use  of  reaction  ( 3 )  to  generate  the 
O  ( ’/*)  atom  quenchers. 1 '  As  this  reaction  takes  place  direct¬ 
ly  in  the  main  flow,  recovery  of  ka  requires  that  the  rate 
coefficient  for  the  reaction 

*NO 

N(*Z>)  +  NO  —  N2  +  O  (4) 

be  known  accurately.  Reported  values'7, 111  for  this  rate  coef¬ 
ficient  range  from  5.9x  10"  "  to  1.8X  10  "’em’ s  thus, 
it  was  necessary  to  carry  out  a  measurement  of  kNO  as  well, 
in  order  to  determine  the  correct  value  of  ka. 

II.  EXPERIMENTAL 
A.  Apparatus 

The  experimental  setup  was  essentially  the  same  as  that 
used  in  previous  rate  measurements  performed  using  the 
FACELIF  reactor14  and  will  therefore  be  described  only 
briefly  here.  A  schematic  is  shown  in  Fig.  I.  The  N(2/)) 
atoms  were  produced  in  a  70  W  microwave  discharge  of  a 
l%-2%  mixture  of  N,  in  helium.  The  N,  and  He  were  puri¬ 
fied  using  gettering  furnaces  containing  Cu  at  600  *C  and  Ti 
at  800  *C,  respectively,  to  remove  O,  impurities  prior  to  en¬ 
tering  the  discharge  tube.  The  NO  (Matheson  98%  purity) 
was  flowed  through  a  purifying  manifold  consisting  of  an 
ascarite  trap  and  a  liquid  nitrogen/methanol  slurry,  and  was 
then  diluted  with  He  to  — 10%  concentration  prior  to  use. 
All  other  quenchers  were  used  without  furthur  purification. 

Metastable  nitrogen  atoms  are  deactivated  at  essentially 
every  wall  collision.* '''  The  N(2Z3)  concentrations  created 
by  microwave  discharges  are  approximately  two  orders  of 
magnitude  lower  than  the  ground  state  nitrogen  atom  con¬ 
centration.  The  combination  of  these  two  factors  made  it 
essential,  under  certain  circumstances,  to  be  able  to  intro¬ 
duce  the  metastable  atoms  into  the  reaction  zone  as  quickly 
as  possible,  while  the  N(:£>)  was  still  in  sufficient  concentra¬ 
tion  to  provide  an  adequate  signal  level.  This  necessitated  a 
fast  pumping  speed  which  was  provided  by  a  Rootes  blower 
( Lcybotd-Heracus  WSI000)  backed  by  a  large  displace¬ 
ment  forepump  ( Heraeus-Englehard  DKI80).  The  maxi- 
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FIC.  I.  Schematic  of  FACELIF  experimental  apparatus.  The  solid  circle 
represents  a  right  angle  connection  to  the  pumping  system. 
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mum  linear  flow  velocity  attainable  using  this  configuration 
was  of  the  order  of  6000  eras'1,  which  corresponds  to  a 
minimum  flow  time  of  approximately  2  ms.  Longer  flow 
times  could  be  achieved  by  reducing  the  pumping  speed  as 
required.  Quenching  gases  were  injected  through  a  loop  in¬ 
jector  whose  position  could  be  adjusted,  over  a  limited  range, 
along  the  length  of  the  flow  tube,  again  altering  the  flow 
time.  Some  comparative  studies  were  performed  using  a  per¬ 
forated  ball  injector  to  probe  any  possible  disparities  in  mix¬ 
ing  properties.  These  indicated  that  the  mixing  characteris¬ 
tics  were  comparable  with  marginally  poorer  mixing  for  the 
ball  injector. 

The  N(2/))  concentration  was  monitored  using  reso¬ 
nance  enhanced  multiphoton  ionization.11'' 14  The  frequency 
doubled  output  of  a  Nd:YAG  pumped  dye  laser  (Quanta- 
Ray  DCR2A,  PDL-I,  WEX-I)  was  used  to  photoionize 
N(2£>),  resonantly  enhanced  by  the  two-photon  transitions 
at  268.93  and  268.98  nm.“  The  current  produced  across  a 
pair  of  Ni  wire  grids,  biased  at  90  V,  was  amplified  using  an 
Ithaco  model  1211  current  preamplifier  and  averaged  on  an 
EG  A  O  model  162  boxcar  averager.  Signal  levels  were  gen¬ 
erally  measured  by  setting  the  laser  on  resonance  at  each 
data  point  and  subtracting  the  off-resonance  baseline  from 
the  averaged  signal.  Some  results  were  also  taken  by  scan¬ 
ning  the  laser  across  the  resonance  lines  with  no  significant 
difference  in  the  results. 

The  method  used  to  generate  the  0(3P)  quenchers  in¬ 
volves  the  titration  of  NO  directly  into  the  flowtube."  The 
NO/He  mixture  is  introduced  into  the  flow,  through  the 
loop  injector,  where  it  reacts  rapidly  with  N(45),21  which  is 
present  in  high  concentrations  in  microwave  discharge  flows 
of  Nj  /He.  Under  ideal  conditions,  the  N(4S)  will  be  in  suffi¬ 
cient  excess  that  all  the  NO  will  be  rapidly  converted  to 
0(3P)  before  making  any  significant  contribution  to  the 
N(2£>)  quenching.  The  diffusion  coefficient  for  O  atoms  in 
He  is  700  cm2  s~  ',22  and  the  surface  recombination  coeffi¬ 
cient  is  on  the  order  of  10  ~  \22  For  the  flowspeeds  and  pres¬ 
sures  used  in  our  2  in.  diameter  flowtube,  we  estimate  less 
than  1%  loss  of  Of2/*)  over  the  reaction  length. 

B.  Experimental  procedure 

Preliminary  experiments  were  carried  out  investigating 
the  variation  in  the  total  N(20)  quenching  rate  over  a  wide 


FIG.  2.  Variation  in  the  effective  rate  coefficient,  k,,  [  Eq.  ( 5)  |  with  reac¬ 
tion  lime.  The  reaction  time  was  varied  by  adjustments  to  the  injector  posi¬ 
tion  and  the  flow  speed.  Consequently,  the  How  conditions  |i.e.,  total  pres¬ 
sure,  |  N(ISj  J  etc.]  varied  with  the  reaction  lime. 


range  of  reaction  times.  These  were  used  to  calculate  an  ef¬ 
fective  rate  coefficient  at  each  reaction  time.  This  is  defined 
as 

*„=  -a(u/rf)(rfln[N(2Z>)]/d[Q]),  (5) 

where  u  is  the  average  flow  speed,  d  is  the  injector  to  detector 
distance,  a  is  the  noniaminar  flow  correction  having  a  nu¬ 
merical  value  of  1.62,*  and  [Q]  is  the  initial  quencher  con¬ 
centration.  The  use  of  Eq.  (3)  is  equivalent  to  assuming  that 
all  the  initially  added  NO  is  converted  instantaneously  to 
OOP)  by  the  N(4S)  without  contributing  to  the  N(2Z>)  de¬ 
cay.  The  results,  shown  in  Fig.  2,  display  a  rapid  decrease  in 
the  effective  rate  coefficient  with  increasing  reaction  time. 
Even  at  ~  50  ms  there  is  some  indication  that  the  rate  coeffi¬ 
cient  is  still  apnroaching  an  asymptotic  value  of  less  than 
1 X  10 “ 12  cm2  s'1.  The  flow  conditions  over  this  range  of 
reaction  times  necessarily  varied  by  a  substantial  amount. 
The  total  pressure  varied  from  0.2  to  4  Torr  in  going  from 
short  to  longer  reaction  times.  There  was  a  corresponding 
change  in  the  N(4S)  concentration,  although  this  was  not 
explicitly  measured  for  all  points  along  the  curve.  The  per¬ 
centage  of  initially  added  NO  remaining  unreacted  at  the 
detection  point  was  also  measured  for  the  same  range  of  re¬ 
action  times.  The  NO  y(0,0)  MPI  signals24-25  were  used  for 
this  purpose.  The  magnitude  of  this  signal  is  directly  propor¬ 
tional  to  the  NO  concentration.  A  blank  titration  was  per¬ 
formed  under  the  same  set  of  experimental  conditions,  ex¬ 
cept  with  the  discharge  turned  off.  The  ratio  of  the  two 
slopes  yields  a  value  of  the  fraction  of  initially  added  NO  that 
remains  unreacted.  The  results  are  shown  in  Fig.  3. 

Comparison  of  Figs.  2  and  3  shows  clearly  that  the  faster 
effective  rate  constants  are  correlated  with  significant  resid¬ 
ual  NO  concentrations.  The  short  reaction  time  measure¬ 
ments  represent  a  region  in  which  reaction  with  NO  domi¬ 
nates  the  N(J£>)  decay.  These  measurements  indicate  an 
extreme  lower  limit  on  Ano.  The  longer  reaction  time  mea¬ 
surements,  to  the  extent  that  there  is  still  some  NO  contribu- 
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FIG.  .V  Variation  in  I  Ik-  fractional  amount  of  initially  added  NO  that  re¬ 
mained  tinreaclcd  at  the  Mri-dclcclion  point.  As  with  the  measure¬ 
ments.  flow  eomlilione  varied  with  reaction  time. 


lion  to  the  effective  decay  rate,  represent  an  extreme  upper 
limit  on  the  0(  7*)  quenching  rate  coefficient.  The  limiting 
values  obtained  Tor  ANO  and  An  obtained  in  this  way.  indi¬ 
cate  a  ratio  Tor  these  two  rate  coefficients  in  excess  of  25  to  I. 
Consequently,  for  short  reaction  times,  where  INOJ  ~  lO], 
the  0(7*)  contribution  is  smaller  by  a  factor  kn /A*,,.  and 
may  be  legitimately  neglected  when  determining  kNO,  Con¬ 
versely,  the  large  value  of  ANO  compared  to  A0  means  that 
even  very  small  amounts  of  NO  may  cause  a  significant  de¬ 
gree  of  quenching. 

It  has  been  reported  that  REMPI  detection  of  N(JD)  at 
269  nm  can  generate  N (*D)  through  the  two  photon  disso¬ 
ciation  of  NO.2*' 27  In  order  to  ensure  that  this  effect  was  not 
compromising  our  results  we  looked  for  the  characteristic 
N (2D)  signal  at  high  NO  flow  with  the  discharge  turned  ofT 
(i.e.,  with  [NO]  significantly  higher  than  normal  operating 
conditions).  No  signal  was  observed  until  [NO]  was  over 
100  times  that  used  during  titrations.  The  laser  power,  mea¬ 
sured  after  exiting  the  flow  tube,  was  <0.5 mJ  pulse"  'com¬ 
pared  to  the  2-5  mJ  pulse  1  reported  by  Jusinski  et  at.1*' 
This  gives  a  production  rate  ratio  of  between  I:l6and  1:100 
for  the  two  systems.  Also.  NO  concentrations  in  the  detec¬ 
tion  region  were  always  <0.5  mTorr.  The  combination  of 
these  two  factors  accounts  for  the  lack  of  any  measurable 
levels  ofNt  '/J)  production  from  NO  dissociation. 

Measurements  of  A  Vl ,  and  A„  both  require  knowledge 
of  [ N(J5)  1  to  determine  how  fast  the  NO  is  depleted  and 
therefore  the  extent  to  which  it  contributes  to  the  N(:D) 
decay.  For  long  reaction  limes  the  standard  NO  titration 
method'’1  may  be  used  to  determine  [N(45)l  accurately. 
This  technique  has  serious  limitations  at  the  shorter  reaction 
limes  to  which  the  NO  rate  coefficient  measurements  were 
constrained.  To  obtain  a  distinct  titration  end  point,  reaction 
( 3 )  must  have  sufficient  time  logo  to  completion.  The  meth¬ 
od  therefore,  is  best  suited  to  conditions,  where  |  N(  ‘ S )  ]  is 
high  and  the  reaction  time  is  long.  This  is  precisely  the  oppo- 
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site  of  the  requirements  for  accurate  kmt  measurements. 
Other  practical  difficulties  were  encountered  in  the  form  of 
greatly  reduced  signal  levels  and  interfering  fluorescence. 
The  fast  flows  and  consequent  low  pressures  (0.2-0. 4  Torr) 
lead  to  significantly  reduced  signal  levels  and  a  subsequent 
reduction  in  signal  to  noise.  In  addition,  stray  light  from  the 
discharge  created  a  background  large  in  comparison  to  the 
target  signal.  Anomalous  effects  were  observed  for  small  NO 
additions  in  the  form  of  an  initial  rise  in  signal  level.  It  is 
possible  that  this  has  its  origins  in  fluorescence  resulting 
from  some  coincident  reaction.  It  has  been  previously  re¬ 
ported2®  that  some  short-lived  source  of  N2  (  B),  other  than 
N(4S)  atom  recombination,  exists  in  discharge  flows  which 
may  interfere  with  titration  measurements. 

Attempts  were  made  to  fit  the  titration  plots  for  short 
reaction  times.  This  introduced  a  further  limitation  beyond 
the  practical  problems  already  mentioned.  Fitting  would  re¬ 
quire  the  use  of  the  rate  coefficient  for  reaction  (3),  which 
has  associated  error  bounds  that  will  propagate  into  the 
[  N(4S)  ]  value.  As  a  result  of  these  limitations,  an  alterna¬ 
tive  method  of  determining  [  N(4S)  ]  was  devised.  Measure¬ 
ments  of  NO  depletion  were  made  using  y  (0,0)  MPI  band 
intensities,  as  described  earlier.  The  high  sensitivity  of  MPI 
to  NO  enables  very  small  concentrations  to  be  detected  and 
it  is  therefore  possible  to  work  with  [N(4S)  ]  >  [NO].  The 
NO  decay  can  then  be  described  by  the  pseudo-first-order 
rate  eqiu  on 

[NO'  —  [NO]0  exp(  —  An  [N(*5)  ]f),  (6) 

where  [K?  _  ]0  is  the  initial  concentration  of  NO  added.  The 
N(4S)  concentration  is  simply  given  by 

[N(4S)J=  —  ln{[NO]/[NO]0}/(ANf)  (7) 

Since  the  ratio  [NO]/[NO]0  may  be  determined  very  pre¬ 
cisely,  the  accuracy  of  this  technique  is  limited  by  the  uncer¬ 
tainty  in  An  . 

For  slow  flows  ( >  20  ms),  it  is  reasonable  to  assume 
that  the  mixing  time  is  short  in  comparison  to  the  total  flow- 
lime.  For  very  fast  flows  (  <  Sms),  this  assumption  may  lead 
to  a  significant  underestimation  of  the  rate  coefficients.  A 
standard  technique  used  to  probe  mixing  effects  is  to  per¬ 
form  rate  measurements  over  a  range  of  injector  positions 
under  constant  flow  conditions.  A  plot  of  the  decay  coeffi¬ 
cient  as  a  function  of  the  flowtime  will  give  a  straight  line 
whose  .r-axis  intercept  gives  a  value  for  the  effective  mixing 
time  and  whose  slope  yields  a  value  for  the  rate  coefficient.  In 
the  case  of  NO.  the  accompanying  reaction  with  N(4S)  de¬ 
pletes  NO  and  thus  compromises  the  technique.  In  order  to 
circumvent  this  problem,  we  performed  the  rate  measure¬ 
ments  substituting  O.  as  a  quencher.  If  it  is  assumed  that  the 
mixing  characteristics  for  O,  and  NO  arc  comparable,  then 
the  effective  mixing  time  obtained  using  O,  may  be  applied 
to  NO.  provided  that  identical  flow  conditions  are  em¬ 
ployed. 

III.  ANALYSIS  OF  RESULTS 

Measurements  of  both  AN(I  and  A„  involve  injecting  NO 
molecules  directly  into  the  main  discharge  flow,  where  they 
react  rapidly  with  N(JS)  atoms  to  generate  0(’/>).  Wlnit- 
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ever  the  flow  conditions,  the  NO  will  be  present  for  some 
finite  time  in  the  tube.  The  N(:X>)  atoms  may  be  quenched 
by  OCP)  and  NO,  both  of  which  will  be  present  in  varying 
concentrations  along  the  length  of  the  Aowtube.  From  reac¬ 
tions  (2)  and  (3),  the  rate  equation  describing  the  decay  of 
N  CD)  is 

-<*[N(2D)1  =  A.  [Nf^KNOj 

dt 

+  *ol  N(lD)][01 

+  Jkw.[N{,2»],  (8) 

where  kw  is  the  first  order  rate  coefficient  for  wall  deactiva¬ 
tion.  Integration  of  this  equation  over  the  reaction  time  gives 

-  ln{[N(J£»l/[N(l2»U  =  kH o  £  [NO)rff 


where  (N(22>)  ]0  is  the  initial  concentration  of  N(JZ3).  The 
concentration  of  O OP)  at  any  given  time  is  simply  equal  to 
[  NO  J , ,  the  amount  of  NO  reacted,  if  there  is  no  other  signifi¬ 
cant  loss  of  OOP). 


The  rate  of  change  of  NO  with  time  is  given  by 

=  k»  [  NO  J  [N(4S)  ]  +  [  NO]  [  N*  ]. 

dt 

(10) 


where  [ N* )  is  the  metastable  nitrogen  atom  concentration. 
Since  [N(4J)J  >  |N*],  the  second  term  in  Eq.  (10)  maybe 
neglected,  and  the  NO  concentration  is  given  by  Eq.  (6). 
Substituting  this  into  Eq.  (9)  gives 


—  In 


[N)(2Z» 

[N(JZ»]0 


=  *olNO]0/  + 


(^NO  ~  ) 

*n[N(4S)] 


x{!  -  exp(  -  kN  [N(4S)]r)} 


X  [NO]„.  (ID 

For  the  situation  where  [  N(4S)  ]  cannot  be  considered  to  be 
in  large  excess,  a  slightly  more  complex  expression  results. 
The  NO  concentration,  in  this  instance,  is  expressed  as 

■  ~-£!.n9L  =  a-  {[NO)0  -  [NO],} 

dt 

x{[N(4S)]0-[NO],}.  (12) 

where  [N^S)],,  is  the  initial  ground  state  nitrogen  atom 
concentration.  This  is  a  second  order  rate  equation  which 
when  integrated  gives” 


1  [[N(4S)]  (tNOl0-(NOl,)  \  # 

([NO]0  — [N(4S)])  "l  [NO]0  <[N(4S)]-[NO],)j 
Rearranging  this  expression  for  [NO]  gives 


[NO],  = 


[N(45)j  -  ([NOUN(45)]/[NO]u)exp{  -  ([NO],, -  [N(4S)])*Nr} 
1  -  ( [N(45) )/[  NO]0)exp{  -  ( [  NO]0  -  [  N(45)  ]  )*N t } 


(13) 


(14) 


The  N(2X>)  concentration  at  any  point  along  the  flowtube  is 
given  by 

xj£  [NO]rff  + const.  (15) 

(NO)df  must  be  evaluated  by  numerical  integration 

over  the  range  /  =  0  to  t^ ,  where  t^,  is  the  reaction  time 
corrected  for  mixing  effects. 

IV.  RESULTS  AND  DISCUSSION 

An  estimate  of  the  effective  mixing  time  for  the  NO  ki¬ 
netic  runs  was  obtained  from  the  O,  rate  measurements.  The 
results  are  shown  in  Fig.  4.  The  total  pressure  in  the  flowtube 
was  —0.4  Torr  with  a  linear  flow  velocity  of  —  6000  cm  s 
A  weighted  linear  regression  lit  to  the  data  yields  a  value  of 
(0.48  ±  0. 19 )  ms  for  the  effective  mixing  time.  The  slope  of 
the  line  in  Fig.  4  gives  a  value  of  the  N(2Z7)  +  O,  rate  coeffi¬ 
cient  of  ( 5.89  ±  0.40)  X 10  cm*  s  \  in  good  agreement 
with  previous  determinations.  ,‘l° 


0  I  i  J  4  S 

T  I  out i  i  me  xms 


FIG.  4-  Variation  in  the  decay  coefficient  (  -  </ln[N('/>)  I  /d  ]0;  ] )  with 
reaction  lime  for  N(  '/»  quenching  by  O.  Flowtube  pre»ure  was  0  46 Torr 
at  a  linear  flow  velocity  of 6700  cm  s  '.I  lie  error  lun  rcpi  cvtii  the  MatiMt- 
cal  limits  to  the  decay  coetTictciu  obtained  from  a  weighted  linear  regression 
til  to  the  individual  data  sets. 
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FIG.).  Decay  of  iheNC'O)  MP1  signal  asa  function  of/;  |NO)rfr,  which  wat  calculated  by  numerical  integration  of  Eq.  (14).  t  m  3.90  ms  (corrected  for 
0.48  ms  mixing  time).  ks  m  J.4x  10  "  cm’  *  and  IN(*S)]  *  0.4.1  mTorr. 


The  NO  experiments  were  carried  out  using  identical 
flow  conditions  to  those  for  the  O,  runs.  The  N  (4S)  concen¬ 
trations  were  varied  between  0.25  and  0.45  mTorr.  For  these 
conditions  f'„  [0]t//~.fJlNO]<//,  but  as  Amo>Ao.  the 
0('/*)  atom  quenching  may  be  neglected.  From  Eq.  (9),  a 
plot  of  the  N( ’/»  signal  decay  vs  [NOJrf/  will  give  a 
straight  line  whose  slope  gives  a  value  for  k directly.  Such 
a  plot  is  shown  in  Fig  5  for ;  =  3.90  ms  and  [  N(45)  J  =  0.43 
mTorr. 

There  arc  a  number  of  sources  of  error  that  must  be 
considered  in  calculating  A***-  Uncertainties  in  signal  levels 
are  shown  as  error  hars  in  Fig.  5.  The  dynamic  range  of  the 
N(:D)  decay  which  could  be  followed  was  limited  by  the 
appearance  of  the  NO  j'(0.4)  band  at  higher  NO  concentra¬ 
tions.  which  obscured  |IkN(jD)  signal  to  some  extent.  This 
baseline  variation  is  reflected  in  the  larger  error  bars  for  high 
[NO]  additions.  The  vibrationally  excited  NO  population  is 
presumably  the  result  of  the  reaction  *" 

N,(/f)  +  NO(.V)--N,(JO  +  NOM)  (16) 

followed  by  j'-band  omission,  although  this  has  not  been  con¬ 
firmed  experimentally. 

The  value  used  for  AN  in  the  calculations  is  taken  from 
the  literature'1  and  has  an  associated  uncertainty  of 
±  26%.  This  value  is  used  in  calculating  f  N(4S) )  and  again 


in  the  integration  of  /„  [  NO)  dr.  This  error  manifests  itself  as 
a  shift  in  the  x-axis  position  of  the  data  points.  The  extreme 
values  for  kH,  [N(4S)  ],  and  were  used  to  calculate  up¬ 

per  and  lower  bounds  for  A*,,.  These  bounds  are  shown  in 
Fig.  6  Tor  the  data  shown  in  Fig.  5.  The  upper  and  lower 
bounds,  together  with  rate  coefficient  obtained  using  the 
mean  values  of  AN ,  [  N(4S)  ]  and  .  are  given  in  Table  I  for 
all  kinetic  runs.  The  error  limits  on  each  value  are  the  statis¬ 
tical  errors  (  ±  I  o)  obtained  from  a  weighted  linear  regres¬ 
sion  fit  to  the  data  points  in  each  limit.  These  arc  small  in 
comparison  to  the  error  bounds  indicating  that  the  domi¬ 
nant  source  of  error  arises  from  the  uncertainty  in  the  values 
of  kK  and  the  mixing  time.  The  final  values  are  obtained  by 
Inking  the  weighted  mean  of  all  kinetic  runs.  This  procedure 
yields  a  result  of  (6.70.  +  1.36,  —  1.21  )X  10  "  cm's  1 
for  kKn . 

The  ideal  conditions  in  which  to  measure  Ay,  are  at  long 
reaction  limes  with  [NO|  <<  [N(45)  J.  If  these  conditions 
can  be  fulfilled  then  Eq.  (II)  may  be  used  to  evaluate  A„. 
Given  the  great  disparity  between  ANO  and  k„,  these  condi¬ 
tions  arc  extremely  difficult  to  realize  in  practice.  At  the 
slowest  flows,  the  use  of  high  flowtubc  pressures  was  not 
sufficient  to  offset  the  loss  of  signal  due  to  increased  wall 
deactivation,  thus  constraining  the  experimental  range  of  i. 
The  maximum  [N(4S)  ]  values  obtained,  even  under  these 
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FIO.  4  Error  bound*  on  the  decay  of  (he  N('fl)  MPI  signal  as  a  function  of  /£(NO)rfi  for  (he  data  shown  in  Fig.  S.  Lower  trace:  >  =  3.71  ms, 
-  4.3X  10  '  “  cm1  s  ~ '  and  1N(\S)  |  .  0.52  mTorr.  Upper  (race:  l  -  4.09  ms.  *„  =  2.5  X  10  "  cm*  s  - '  and  |  N(4£) )  .  0.34  mTorr. 


extreme  conditions,  were  1.6XI014  cm" 1  ( ~5  m  Torr). 
The  low  value  of  k0  required  NO  additions  of  up  to  3.5 
mTorr  to  give  (N(aA))/[N(,/>))0  <0.05  and  thus  the 
(N(4S) ] >  [NO]  condition  could  not  be  met.  Some  at¬ 
tempts  were  made  to  enhance  N-atom  dissociation  by  add- 


TABLE I.  Kinetic  data  for  determination  of  N('f>)  +  NO-N,  +  O  rale 


• 

/,  ms 

Lower  bound 

molecule' 1  s ' 
Mid  range 

•  xio  " 

Upper  bound 

I.S2 

6.11  ±0.18 

6.63  ±  0.2S 

7.72  ±  0.22 
8.52  ±0.36 

8.79  ±0.26 
9.73  ±0.41 

2.625 

7.3S  ±  0.31 

7.53  ±  0.30 

8.77  ±0.37 
1.88  ±  0.35 

10.30  ±0.43 
10.60  ±0.42 

3.175 

5.81  ±0.22 

6.64  ±0.40 

7.13  ±0.27 
8.24  ±  0.50 

8.89  ±0.34 
10.40  ±  0.63 

3.315 

6. 56  ±0.30 

6.56  ±0.26 

7.84  ±0.36 
7.76  ±0.31 

9.47  ±0.43 
9.28  ±  0.37 

3.165 

4.30  ±0.11 

4.61  ±  0.24 

3.13  ±0.22 
5.31  ±0.28 

6.06  ±  0.23 
6.38  ±0.33 

3.895 

Weighted 

4.61  ±0.13 

5.00  ±  dll 

5.67  ±0.16 
6.13  ±0.22 

7.13  ±0.20 
7.63  ±  0.27 

SMM 

5.49  ±0.07 

6.70  ±  0.01 

8.06  ±0.10 

ing  SF*  to  the  discharge  flow.11  This  did  result  in  a  signifi¬ 
cant  increase  in  [N(45)]  (by  a  factor  of  up  to  2.5)  and  also 
in  the  N(  JZ>)  signal  level.  Unfortunately,  the  use  of  SF6  also 
resulted  in  the  appearance  of  a  large  molecular  MPI  spec¬ 
trum  in  the  same  frequency  range  which  also  varied  with 
quencher  concentration  and  therefore  compromised  the  re¬ 
sults.  It  is  possible  that  such  interferences  would  not  occur 
for  other  N(2Z>)  detection  methods.  If  this  were  true,  then 
this  method  could  be  used  to  enhance  both  N(4$)  and 
N(2Z)),  assuming  that  the  species  generated  by  discharging 
SFt  do  not  participate  in  the  chemistry  in  any  way. 

Equation  (15)  was  used  in  calculating  ka.  The  pre¬ 
viously  determined  value  of  k *<,  was  used  and  the  NO  con¬ 
tribution  to  the  decay  subtracted  at  each  point.  Figures  7  and 
8  show  the  efTect  of  applying  this  correction  to  the  raw  data 
for  two  separate  kinetic  runs.  The  relative  contribution  from 
NO  increases  with  increased  quencher  addition  and  shorter 
reaction  time.  Even  for  the  highest  [N(4£)]  values  and 
longest  flowtimes,  NO  accounted  for  —30%  of  the  N(  JZ>) 
decay. 

The  N (*S)  concentration  can  be  measured  by  the  more 
usual  titration  technique  for  these  slow  flow  conditions. 
[N(4S)]  was  measured  before  and  after  each  kinetic  run 
with  results  agreeing  to  within  2%  in  all  instances.  We  con¬ 
servatively  estimate  an  accuracy  of  ±  5%  for  the  NO  titra¬ 
tion  technique.  The  mixing  time  was  not  explicitly  measured 
for  these  longer  reaction  times  as  this  can  reasonably  be  ex- 
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FIG.  7.  Decay  of  the  N('0)  MFI  signal  as  a  function  of/;  (Oldf.  The  bottom  trace  was  calculated  assuming  instantaneous  conversion  of|NO]  lo(01  (ii, 
r.(0|4r-  INOl.')  Md  therefore  no  correction  for  depletion  of  N(’D)  by  NO  h  necessary.  A  linear  regression  «t  to  these  data  gives  a  value  of 
ho  —  2.0)  X 10' 11  an’ t' '.  The  upper  trace  was  calculated  using  Eq.  (15)  with  *M  -  J.4x  10'"  an’ s*'.**,  -6.7X10*”  cm's-',  |N(*5))  -  3.21 
mTorr  and  r»  23.0  ms.  This  gives  a  value  of  «■  7.40X  H>  *  “as1 «  The  difference  in  fhetwpplotsshows  the  ndaiive  magnitude  of  the  NOccntiftmioa 
to  the  N('D)  signal  decay. 


/  (03. dt  xl0l<cm  ’a 

FIO.  ».  tame  as  Fig.  7  for  r  ■»»  41.4  mv  The  top  I  race  was  calculated  asing  —  3.4  x  10  "em’s  fc.Tx  10  "cm's  1  and  |N(*S)  I  -■*  4.** 

mTorr. 
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}  CO)  .dt  x  10'*cm  *s 


FIG.  9.  Upper  and  lower  bounds  to  N('O)  decay  as  a  function  of  f‘, (O)rff  for  Ihe  data  shown  in  Fig.  7.  Lower  trace:  ks  =  4.Jx  10  "  cut's 
A*,,  *  8.06  x  10  "  cm's" 1  and  |N(*S)|  »  3.44  mTon.  Upper  trace:  kK  =  2.5x  10  "cm's  5.49  x  10  "cm's  1  and  fN('.?)|  =  3.12 

mTorr 


peeled  to  be  short  compared  to  the  flowtime.  Given  the  long 
flowtime  and  large  N(4.S)  concentration, 
io  [NO]u/  and  the  correction  for  the  NO  contribu¬ 

tion  to  the  decay  involves  primarily  a  y-axis  adjustment.  The 
magnitude  of  this  adjustment  is  sensitive  to  the  exact  values 
used  for  kN  and  Ar^o-  In  the  interests  of  self-consistency 
within  the  analysis,  when  one  extreme  value  of  kH  was  taken 
in'  the  error  calculation,  the  value  of  k^,  calculated  using 
that  same  value  was  used.  Thus  the  kN  and  k^  errors  are 
correlated  and,  in  this  instance,  tend  to  cancel  one  another. 
The  measurement  of  kit  is  therefore  less  sensitive  to  the  un¬ 
certainty  in  kN  than  is  ftNO,  The  error  bounds  for  k„  were 
calculated  using  the  extreme  values  of  (N(45)|,  kN,  and 
kNu.  The  results  are  shown  in  Figs.  9  and  10  for  the  data  in 
Figs.  7  and  8,  respectively.  The  results  of  twelve  separate 
kinetic  runs  are  given  in  Table  II.  As  before,  the  errors 
quoted  for  each  value  are  the  statistical  values  in  each  limit 
and  are  small  compared  to  the  uncertainties  arising  from  the 
correlated  systematic  errors.  The  weighted  averages  of  these 
data  give  a  value  for  kQ  of  (6.93  4-  0.71,  —  1.08) X 10"  13 
cm1  s  " 1 .  The  absence  of  any  systematic  variation  in  k„  with  t 
confirms  that  there  is  no  significant  loss  of  O (3P)  along  the 
flow  tube.  The  02  formed  would  react  rapidly  with  N(20), 
giving  an  apparent  increase  in  k0  at  longer  reaction  times. 


V.  CONCLUSIONS 

The  previously  reported  experimental  and  modeling 
values  of  k0  are  listed  in  Table  III.  The  value  of  ka  reported 
here  is  in  good  agreement  with  the  most  recent  experimental 
measurement  made  by  Piper,32  who  utilized  Ol  'P)  atom 
generation  and  N(2£)  detection  techniques  differing  from 
those  employed  here.  It  is  also  in  fair  agreement  with  the 
earlier  results  of  Davenport  et  al.}i  and  lanuzzi  et  al ."  The 
major  conflict  is  with  the  result  of  Jusinski  et  at.,"  who  re¬ 
cently  rekindled  the  ongoing  controversy  regarding  the  dis¬ 
parity  between  experimental  measurements  and  the  pre¬ 
ferred  modeling  values  by  suggesting  that  the  true  rate  was  of 
Ihe  order  of  2.1  X  10'"  cm’s' '.  The  experiment  by  Ju¬ 
sinski  et  al.  and  the  studies  reported  here  were  performed 
using  an  essentially  identical  experimental  setup  and  meth¬ 
odology.  The  origin  of  the  disparity  in  the  results  is  not  im¬ 
mediately  apparent  although  it  is  interesting  to  note  that  our 
preliminary  investigations  (as  illustrated  in  Fig.  2)  also  yield 
a  value  of  ~  2  X 10 '  1 1  cm' s ' 1  at  /  =  7. 5  ms.  Our  results 
and  calculations  indicate  that,  for  the  conditions  reported  in 
Ref.  11  (lN(4S)]  =  lOmTorr.f  =  7.5  ms),  an  effective  rate 
constant  (which  is  what  was  actually  reported)  of 
~2x  10' 12  cm’s'1  should  be  observed.  The  order  of  mag¬ 
nitude  difference  cannot  be  accounted  for  by  any  reasonable 
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JcO].dt  x  I0'*cm  *s 


FIG.  10.  Upper  and  lower  bounds  for  N(J2>1  decay  as  a  function  of  p„\0\dl  for  the  data  shown  in  Fig.  8.  Lower  trace:  ku  — 43  x  10  "  cm's'’. 
*,,,  -8.06x10  "  cm’s  *  and  IN(*S)|  -4.49  mTorr.  Upper  trace:  A*  -2.5x10'  "  cm’s  =  5.49x10  "  cm’s  and  |N(*S)|  -40| 

mTorr. 


adjustment  for  poor  mixing.  It  is  possible  that  the  faster  rate 
observed  is  the  result  of  the  quenching  of  some  precursor 
which  is  the  dominant  N  ( 20)  source  in  their  higher  pressure 
system  ( ~30  Torr).  The  nonlinearity  observed  in  some  of 


TABLE  II.  Kinetic  data  for  determination  of  N('l>)  +  Ot  'P)  rate  coeffi¬ 
cient. 


cm ’ molecule'1  *  ’x  10 

-  II 

/,  ms 

Lower  bound* 

Mid  range’’ 

Upper  bound* 

23 

4.74  ±  0.37 

7.60  ±0.37 

9.17  ±0.37 

3.37  ±0.50 

6.38  ±  0.49 

8.03  ±  0.49 

317 

5.33  ±  0.23 

6.81  ±  0.23 

7.57  ±0.23 

5.01  ±0.31 

6.78  ±0.31 

7.64  ±0.31 

361 

7.12  ±0.26 

8.51  ±  0.26 

9.27  ±  0.26 

6.01  ±  0.24 

7  44  ±  0.24 

8.27  ±  0.24 

41.4 

6.62  ±0.21 

7.58  ±0.21 

8.14  +  0.21 

6.01  ±  0.26 

6  99  ±0.26 

7.53  ±0.26 

44.8 

5.46  ±  0.28 

6  21  ±  0.27 

6.63  ±  0.27 

5.32  ±  0.28 

6.09  ±  0.28 

6.52  ±  0.28  . 

55 

6.19  ±0.17 

6.65  +  0.17 

6.90  ±  0.t7 

5.54  ±  0.23 

6.05  ±  0.23 

6.33  ±  0.23 

Weighted 

mean 

5.85  ±  0.08 

6.93  ±  0.08 

7.64  ±0.08- 

*U«og  kn 

-2,5X10  "cm 

'*  -5.49X10  " 

cm  *  %  ’  V 

*  Uiing  kM 

-3.4x10  "cm 

V  '.  *w>  -  6.70X  10  " 

cm'*  '. 

■  U*ing  ** 

-  4.3 X  10  "cm1 

i  '.  -=  8  06x10  " 

cm*  *  V 

their  plots  would  tend  to  support  this  hypothesis. 

The  method  of  generating  O CP)  atoms  necessarily 
means  that  NO  molecules  will  be  present  in  the  flow  for  some 
finite  time.  As  k ^  turns  out  to  be  nearly  2  orders  of  magni¬ 
tude  larger  than  ko.  this  will  always  be  significant  although, 
by  carefully  controlling  the  experimental  conditions,  the  NO 
contribution  can  be  kept  to  a  minimum.  It  is  essential  that 
the  flow  conditions  be  precisely  characterized  for  each  kinet¬ 
ic  run,  and  that  the  analysis  carried  out  with  extreme  care,  if 
meaningful  results  are  to  be  recovered.  This  also  applies  to 
the  NO  results.  In  this  instance,  the  requirement  that  the 
measurements  be  carried  out  at  very  short  reaction  times 
brings  us  to  the  limits  of  flow  tube  technology.  The  mixing 
time,  using  standard  injectors  and  flows,  becomes  a  signifi¬ 
cant  fraction  of  the  total  flow  time  and  must  be  carefully 
considered  to  avoid  reporting  artificially  low  values.  The  ac¬ 
curacy  of  the  reported  value  of  kNO  is  primarily  dictated  by 
the  confidence  limits  associated  with  kN .  Should  a  more  pre¬ 
cise  value  of  kN  become  available,  the  confidence  limits  on 
kN 0  may  be  reduced  accordingly. 

Experimental  measurements  of  *o  encounter  a  number 
of  difficulties,  the  most  immediate  of  which  are,  the  genera¬ 
tion  of  a  pure  source  of  0(  'P)  atoms,  their  transportation  to 
the  reaction  zone  and  the  subsequent  efficient  mixing  of  the 
reactants.  The  method  of  titrating  NO  directly  into  the  flow 
tube  ingeniously  solves  the  transportation  and  mixing  diffi¬ 
culties  but  introduces  an  additional  complication  as  the 
0(V*)  precursor  itself  is  a  highly  reactive  impurity.  With 
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TABLE  III.  Experimental  aad  a  representative  sample  of  mode  tin*  determinations  of  *„ . 


*.,cmJs-' 

Source 

Reference 

Experimental 

( 1.8*  0.4)  x  10  12 

Davenport  ti  at.  ( 1976) 

33 

<1.8x10- 12 

Ian  uni  and  Kaufman  ( 1980) 

9 

(2.1±Q.8)X10  " 

Jusinski  rial.  (1988) 

II 

<  1.06  ±  0.26)  x  10  12 

Piper  (1988) 

32 

(6.9*  IO)xlO  " 

This  work 

Modeling 

IXlO12 

Strobel  et  el.  (1970) 

2 

<10  M 

Oran  era/  (197$) 

12 

4X10  12 

Frederick  and  Rusch  ( 1977) 

13 

IX  10 

Cravens  rt  at.  ( 1979) 

14 

4  X 10" 12 

Rusch  and  Gerard  ( 1980) 

1$ 

6X  10" 12 

Richards  rial.  (1981) 

36 

-IXlO'1 2 * 4 

_ 

Fesen  Hat.  (1989) 

3$ 

precise  characterization  of  flow  conditions  and  careful  anal¬ 
ysis  of  results,  an  accurate  value  for  kQ  can  be  recovered. 

Recent  models  for  atmospheric  NO  concentrations 
have  tended  to  favor  larger  values  of  k0 ,  and  it  would  appear 
that  experiment  and  theory  have  finally  reached  a  consensus 
regarding  the  preferred  value  of  this  rate  coefficient.  Bates14 
has  shown  that  the  fast  rate  coefficient  reported  by  Jusinski 
et  a/."  is  incompatible  with  observed  data  on  Ol  red  line 
emission  at  630  nm.  Also,  the  most  recent  modeling  calcula¬ 
tion  by  Fesen”  successfully  accounts  for  observed  NO  dis¬ 
tributions  using  a  value  of  kQ  *  IX  10~,J  cm’s"',  which 
coincides  with  the  most  recent  experimental  determinations. 
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DETECTION  OF  Nf*S)  BY  RESONANTLY  ENHANCED  MULTIPHOTON  IONIZATION  SPECTROSCOPY 


(Received  l  December  1989) 


Abstract  •  We  report  the  first  instance  of  detection  of  N(4S)  atoms  by  Resonantly  Enhanced 
Multipboton  Ionization  (REMPI).  The  magnitude  of  the  ion  signal.  after  correction  for  laser  power 
dependence,  is  directly  proportional  to  the  N?S)  concentration  and.  once  calibrated,  is  a  direct  and  highly 
sensitive  method  for  N(*S)  concentration  measurements  in  fast-flow  kinetics  experiments. 


INTRODUCTION 

In  experiments  utilizing  active  nitrogen 
flows,  it  is  generally  necessary  to  know  the 
ground  state  nitrogen  atom  concentration.  This 


value  is  usually  delenninad  using  a  standard  NO 
titration  technique  [1].  This  method  makes  use  of 

the  reactions 

Nf4S)  +  NO-»N2+0(3P) 

(1) 

and 

NO+O+M  -» NOz’  +  M 

(2) 

Reaction  (1)  is  rapid  [2]  and  serves  to  deplete 
{N(4S)J  stoichiometricslly.  Monitoring  the  Nj 
first  positive  emission  around  580nm  gives  a 
fluorescence  signal  proportional  to  (N^S)]2  for 
[NOJ<(N(4S)].  Beyond  the  end  point  NO2' 
emission  following  reaction  (2)  gives  a  signal 
proportional  to  NO  in  the  same  spectral  range. 
This  technique  works  well  provided  that  reaction 
(1)  has  sufficient  time  to  go  to  completion  and 
N(4S)  is  present  in  sufficient  quantity.  If  this  is 
not  the  case  then  both  the  N2  first  positive  and 
NOj*  emission  will  contribute  to  the  signal 
level. 

Many  reactions  involving  m  etas  table 
atoms  and  molecules  require  fast  flow  conditions 
due  to  the  rapid  decay  of  the  metastable 
concentration.  The  reactions 


N(JD)+N0-tN2+O 

(3) 

Nf^+Oz-tNO+O 

(4) 

are  examples  of  this  [3-3].  In  order  to  calculate 
the  rate  coefficient  of  reaction  (3),  (N(4S)J  must 
be  known  in  order  to  determine  the  rate  of 
depletion  of  NO  via  reaction  (l).  Under  certain 
conditions,  the  flow  time  and  N(4S)  concentration 
may  be  such  that  reaction  (1)  does  not  go  to 
completion  and  the  standard  NO  titration 
technique  will  no  longer  give  accurate  values  for 
(N(4S)].  The  use  of  Resonantly  Enhanced 
Multiphoton  Ionization  (REMPI)  to  measure 
relative  N(^D)  and  N(^P)  concentrations  is  well 
established  (3,6).  Here  we  report,  to  our 
knowledge,  the  first  instance  of  REMPI  detection 


of  N(*S)  atoms.  The  magnitude  of  the  ion  signal 
is  directly  proportional  to  the  N(4S  ) 
concentration  aid,  once  calibrated,  may  be  used 
10  measure  Nf*S)  absolutely. 

EXPERIMENTAL 

These  investigations  were  carried  out 
using  the  FACEUP  afterglow  reactor  which  has 
been  described  in  detail  elsewhere  [4.7],  A 
schematic  of  the  experimental  apparatus  is  shown 
in  Fig.  1.  The  N(*S)  atoms  are  generated  in  a  flow 
of  1-2%  N2  in  He,  using  a  microwave  discharge. 
The  current  created  by  the  photoions,  across  a 
pair  of  Ni  wire  grids,  biased  at  90V.  is  amplified 
using  an  Ithaco  model  1211  current  amplifier  and 
detected  on  an  EG  &  G  model  162  Boxcar 
Averager.  The  laser  power  is  also  monitored 
using  a  Molectron  Pyroelectric  joulemeter  (13- 
09).  The  fluorescence  measurements  are  made 
using  a  photomultiplier  tube  through  a  380  ±  S 
ran  filter.  The  NO  used  for  these  titrations  is  a  5% 
mixture  of  NO  in  He  and  is  added  through  a  loop 
injector  upstream  of  the  detection  region. 
Typical  flow  times  were  25-30ms  with  total  tube 
pressures  of  2.5-3.5  Tore.  These  slow  flow/high 
pressure  conditions  are  best  suited  to  making  NO 
titration  measurements  of  (N(4S)]  giving 
comparatively  intense  signals  and  allowing 
reaction  (1)  ample  lime  to  go  to  completion.  The 
N(4S)  concentration  could  r  t  varied  either  by 
adjusting  the  N2  to  He  ratio  tihis  was  typically  1  - 
3%  Nj  in  He)  or  more  conveniently  by  altering 
the  microwave  power. 

The  transition  N(2p*3p  4Dj°) «-  N(2p^ 
4S]/2°)  may  be  resonantly  excited  by  two  21 1  nm 
photons  (8),  and  the  atom  then  ionized  by  a  third 
photon  of  the  same  energy.  Laser  radiation  at 
211  nm  may  be  generated  by  stimulated  Raman 
frequency  shifting  of  the  doubled  output  of  a 
Nd:YAG  pumped  dye  laser  [9]  but  this  generally 
gives  pulse  energies  of  less  than  100  mJ. 
Significantly  higher  pulse  energies  are  desirable 
for  REMPI  detection.  Here  we  have  used  a  BBO 
crystal  to  double  the  output  from  a  Nd:YAG 
pumped  dye  laser  (Quanta-Ray  DGR2A.  PDL-1, 
INRAD  Autotracker  II)  directly.  Stilbenc  420 
laser  dye  gave  pulse  energies  of  approximately 
0.3  mJ,  measured  after  exiting  the  flow  tube. 
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Thu  power  level  decreased  rapidly  in  time  with 
the  dye  having  a  useful  lifetime  of  about  20 
minutes.  The  decrease  proved  advantageous  when 
taking  measurements  of  the  power  dependence  of 
the  ion  signals,  enabling  a  large  number  of  data 
points  to  be  taken  without  the  extensive  use  of 
neutral  density  filters  to  attenuate  the  beam. 
When  making  quantitative  measurements  of 
(N(^S)l ,  however,  it  is  desirable  to  keep  the  laser 
power  as  constant  as  possible  to  minimize 
corrections  necessary  for  power  normalization  of 
km  signals.  We  found  that  the  use  of  -  four  times 
larger  dye  reservoirs  greatly  enhanced  the  useful 
life  of  the  dye  and  also  minimized  power 
variations  during  the  course  of  data  acquisition. 

ABSOLUTE  (N(4S)1  MEASUREMENTS 

The  N(4Dj°)  state  is  split  into  four 
levels.  1-1/2.  3/2.  3/2.  and  7/2.  at  94770.83. 
94793.46.  94830.86.  and  94881.79  cm1, 
respectively.  The  relative  intensities  are  simply 
given  by  the  upper  state  degeneracy  (2J+1)  and 
the  observed  spectrum  will  be  a  quartet  with 
intensity  ratios  of  1:2:3:4  (J=»l/2:3/2:5/2:7/2). 
The  observed  spectrum  is  shown  in  Fig.  2.  The 
most  intense  line  was  used  throughout  for  all 
quantitative  measurements  of  (N^S)].  To 
normalize  the  signal  levels  accurately  for 
variations  in  laser  power,  it  was  necessary  to 
characterize  the  power  dependence  of  the  ion 
signal  precisely  prior  to  any  data  acquisition. 
The  ion  signal  is  proportional  to  some  power  N 
of  the  laser  intensity  /.  and  therefore  a  plot  of  In 
(ion  signal)  vs.  In  /  will  have  a  slope  equal  to  N. 
The  results  of  doing  this  are  shown  in  Fig.  3.  The 
slope  varies  between  2.8  and  1.3  going  from  low 
to  higher  power.  This  indicates  that  the  two 
photon  transition  is  becoming  saturated  at  the 
higher  powers.  A  similar  power  dependence  was 
observed  for  (2+1)  photon  ionization  of  N(*D)  at 
269nm  [10],  although  somewhat  different 
conclusions  were  drawn  in  that  work  regarding 
the  relative  saturation  of  the  excitation  and 
ionization  steps.  The  two  photon  excitation 
cross-section  for  this  transition  has  been 
calculated  [9]  and  therefore  the  transition  rate 
may  be  calculated.  The  fit  of  the  data  in  Fig.  3  is 
obtained  using  a  simple  rate  equation  model  for  a 
three  level  system  (11].  Strictly  speaking, 
correct  normalization  of  signals  for  power 
variations  requires  that  signal  and  laser  power  be 
recorded  for  each  pulse  (12).  Here  we  have 
recorded  the  average  ion  signal  and  laser  power 
over  a  large  number  of  pulses  with  the  beam 
tightly  focused  between  the  detector  grids.  The 
resulting  ion  signal  therefore  consists  of 
contributions  over  a  range  of  saturation 
conditions.  It  is  thus  inappropriate  to  comment 
on  the  magnitudes  of  the  experimentally 
observed  rates  for  either  the  ionization  or  'be  two 
photon  excitation.  For  the  purposes  of  this  study 
it  is  only  necessary  to  be  able  to  understand  the 
bulk  behavior  of  the  system.  Fig.  3  indicates 
that  for  small  variations  in  average  laser  power 
(±3%)  it  is  a  reasonable  approximation  to  use  the 


same  laser  power  dependence  factor  when 
normalizing  the  ion  signals. 

For  conditions  where  reaction  (1)  has 
time  to  go  to  completion.  (N(4S)]  may  be 
determined  by  observing  the  reduction  in  the  ion 
signal  as  NO  is  added.  The  x-axis  intercept  on  a 
plot  of  the  ion  signal  vs.  [NO]  gives  the  value  of 
[N(^S)1.  The  slope  of  this  line  serves  to  calibrate 
the  ion  signal  level  in  terms  of  absolute  (N(4S)| 
concentrations.  In  order  to  verify  the  accuracy  of 
this  technique  for  (N(4S)]  measurement, 
comparative  measurements  were  taken  using 
REM  PI  and  the  standard  fluorescence  monitoring 
titration  technique  (1].  The  results  are  listed  in 
Table  1.  The  REMPI  results  are  the  x-axis 
intercept  of  a  weighted  linear  regression  Tu  of  NO 
titration  data  sets  with  ±20  error  limits.  The 
standard  titration  measurements  may  be 
considered  accurate  to  ±3%.  The  results  are  in 
good  agreement  indicating  that  the  REMPI 
technique  yields  accurate  values  for  (N^)]. 

DISCUSSION 

REMPI  provides  a  sensitive  method  for 
quantitatively  detecting  ground  state  nitrogen 
atoms.  Number  densities  on  the  order  of  1010 
atoms  cm"3  were  easily  detected  compared  to  a 
limit,  for  our  flow  system,  of  around  2x10' 3 
atoms  cm"3  using  the  more  usual  titration 
technique.  Below  this  level  the  fluorescence 
measurements  encountered  problems  in  the  form 
of  stray  light  from  the  discharge  reaching  the 
detector.  A  mom  serious  problem  can  arise  from 
fluorescence  due  to  sources  other  than  N(*S)  atom 
recombination,  which  contributes  significantly 
to  the  signal  intensity  (4.13].  Unlike  the 
standard  technique,  REMPI  offers  a  state  specific 
method  of  monitoring  (N(4S)]  directly  and  is 
therefore  less  susceptible  to  extraneous  signals. 
The  signal  level  may  be  calibrated  under  slow 
flow/high  (N(4S)]  conditions.  The  magnitude  of 
the  ion  signal  may  then  be  used  to  give  a  direct 
measurement  of  INC'S)]  for  experiments  requiring 
fast  flows  where  (N(4S))  may  be  low;  these 
condition  are  inappropriate  for  standard  titration 
measurements. 
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Table  1.  Comparison  of  (N(4S)]  titrations  by  NO  using  standard  chemiluminescence  detection  and  direct 
MPI  detection  under  flow  conditions.  (N(4S)]  given  in  mTocr. 


Standard  NO  Titration 


MPI -Titration 


4.38 

3.28 

165 

2.02 

1.97 

1.42 


4.48  ±0.18 
3.4010.10 
2J010.14 

1.9610.12 

1.9710.12 
1.5110.07 


» 


FIG.  1.  Schematic  of  FACELIF  experimental  apparatus.  The  solid  circle  represents  a  right  angle 
connection  to  the  pumping  system. 
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I  n  C ION  SIGNAL  (arb.  units)] 


